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Abstract. Using several illustrative examples, the nature of resonance poles and the corres- 
ponding zeroes of the s-wave S matrix is examined for several potentials having an absorp- 
tive pocket followed by a barrier. It is shown that even though the presence of absorption practically 
suppresses the manifestation of resonance in the elastic scattering cross section, the effect of 
the resonances generated by the absorptive pocket is more clearly manifested in the absorption 
cross section provided the barrier width is not too large. We further find that the signature of 
barrier top resonances are also more clearly manifested in the absorption cross section rather than 
in the elastic scattering cross section. These results have been interpreted in terms of complex 
resonance poles and corresponding zeroes of the S matrix. This implies that in complex potential 
scattering like heavy ion collisions, the reaction channel cross section peak is a more reliable 
signature of resonance phenomenon than the variation of the elastic channel cross section with 
energy. 
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1. Introduction 

The analytic properties of the S matrix for real well-behaved potentials, their applications 
in scattering theory and the interpretation of physical states of quantal systems, are well- 
established subjects [1-3]. However, in nuclear scattering, the effective potential between 
the interacting systems is complex and detailed investigations and applications of the 
analytic properties of the S matrix for complex nuclear potentials have been studied only 
to a limited extent [4-7]. Recently, there has been some renewed interest in the nature of 
the singularities of the S matrix for complex potentials and its applications to nuclear 
physics [8-11]. It may be noted that in complex potential scattering, the primary 
observables of interest are the scattering cross section cr s and the absorption cross section 
(Tab- It is desirable to study to what extent both of these quantities can be correlated with 
the analytic features and poles of the S matrix. 



We use the convention 2m = 'h 1 and energy E = k 2 . In numerical calculations, we 
express energy in fm~ 2 unit and length in fm. The lost functions /(fc) are defined as the 
Wronskian W[<f>(k,r} t f(k,r)] = f(k} where <j>(k,r) is the regular solution of the 
modified radial Schrodinger equation with potential U(r}, 

=0 (1) 



having the behaviour $(r) ~ r near the origin and the Jost solutions /(dbfc, r) of eq. (1) 
behave as Lt r ^oo exp(rBfcr)/(fc, r) = 1. Some authors, as a matter of convention use the 
symbol f(k} for f(k) and vice versa (for example, see refs [1] and [8]). 

For real potentials, the Jost functions satisfy the symmetry relations /*(&) =/(&*) 
which means that for real k,f*(k} =f(k}. We denote any typical pole position in the k- 
plane by k P and the corresponding zero by k z . The analytic properties of f(k} ensures 
that if kp = Re(fcp) - i Im(/c/>) is a zero of f(k), then the corresponding zero of f(k} is 
k z -= k* p . This makes the distribution of the poles and corresponding zeroes of the S 
matrix symmetric with respect to the real k axis. For potentials with a pocket, resonance 
states occur inside the pocket which can decay by tunneling through the barrier between 
the outer region and the pocket. For large barrier parameters, the imaginary part of the 
resonance poles become exceedingly small and this can lead to very sharp and narrow 
resonances implying that extremely high resolution is needed to resolve them. However, 
when an absorptive complex potential 

V(r) = U g(r} (2) 

with /o complex having Im(i/o) < and form factor g(r} is used, the symmetry relation 
of the Jost function changes [4] to/*(&, C/o) =/(&*, UQ}. As a result, the zeroes and 
poles of the 5 matrix are no longer symmetrically distributed with respect to the real k 
axis. Hence, we cannot conclude that there is a zero of the S matrix at k k* p 
corresponding to a pole at k kp. Further, resonance states inside the pocket can have 
larger widths of the order of Im (C/o) or larger and hence can give significantly broader 
resonances. It is therefore important to see to what extent such resonances get manifested 
in cr s and <r a b as a function of energy. In the case of complex potential scattering, the 
zeroes of the S matrix have additional significance due to the fact that a^ is governed by 
the factor (1 |S(fc)| 2 ). This means that one can expect a a b to have a maximum when 
\S(k}\ is a minimum. This minimum can be expected to occur in the vicinity of the 
complex zero of the S matrix. Thus peaks in cr ab are related in this sense to the zeroes of 
the S matrix. 

This raises the question as to whether the resonance poles which should get manifested 
in cr s (which gets suppressed due to the large Im (&/>)) can be related to the peaks in cr ab 
which occur due to the zeroes of the S matrix. Such a clear correlation can come only 
when Re (kp} and the corresponding Re (kz} are same or very close to each other. In this 
paper, we examine this problem using several examples and demonstrate that even though 
the peaks in a s get suppressed due to the presence of the complex potential, the presence 
of these states are manifested in <r a b at about the same energies when the imaginary part 



because the real part 01 the resonance poles and the corresponding zeroes are round to be 
very close to each other in the energy plane. This means that it is appropriate to look 
forward to the nature of the resonance states by studying the peaks of the absorption cross 
section as a function of energy. In 2, we describe the details of our calculation and 3 
gives the discussion and conclusion. 

2. Resonance poles and zeroes of S matrix for complex potentials 

In this section we calculate the zeroes and poles of the s-wave S matrix for three different 
absorptive potentials having a barrier in the outer region. Such potentials can generate 
resonance states of comparatively large widths inside the potential pocket. We then 
examine the variation of a s and <r a b as a function of E and interpret them in terms of the 
poles and zeroes of the S matrix. 

2. 1 Rectangular barrier with absorptive pocket 

We first study the poles and zeroes of the S matrix of the potential 

U(r] = -iW, r<a 

= U B: a < r < b (3) 

= 0, r > b 

with 1/5 = 8 fm~ 2 , W = 0.2 fm~ 2 , a = 5 fm, b = 5.5 fm. This gives an absoiptive pocket 
in the region r < a, followed by a barrier of height UB and width (b a). The s-wave 
S matrix is obtained analytically by solving the radial Schrodinger equation which 
is straightforward and we omit the algebraic details. The poles and zeroes of the S-matrix 
are calculated using Newton-Raphson iterative procedure. In figure la, we show the 
variation of <T S and <r ab for s-wave as a function of energy. Curve B indicates cr s obtained 
with W = 0, i.e. curves C and B indicate the scattering cross section in the presence 
and absence of absorption inside the potential pocket. The minima of cr s are very small 
but non-zero. The minima in cr s are due to the variation of the hard core type phase 
shift present in the S matrix and are not intimately correlated with resonance. Figure Ib 
shows the resonance poles and corresponding zeroes in the E-plane. Table 1 lists 
the location of the same poles and zeroes in the fc-plane. The points in the upper (lower) 
half k 2 -plane in figure Ib should be understood to belong to the first (second) Riemann 
sheet. Same observation holds for figure 2b and figure 3b described later. 

It is clear that the introduction of the imaginary part in the potential practically 
eliminates the sharp resonances of curve B from curve C. It should be understood that the 
broad peaks in cr s are not the state resonances but arise as shape resonances. We also note 
that the location of the peaks in curve A are practically at the position of the peaks in 
curve B. The corresponding peaks in curve C are unable to manifest itself due to the large 
width arising from the complex potential. 
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Figure 1. (a) Variation of cr s and cr a b with energy for the rectangular barrier of eq. 
(3) with U B = 8 fm~ 2 , W = 0.2 fm~ 2 , a = 5 fm and b = 5.5 fm. Curve A is the plot of 
cr a b and curves C and B correspond to cr s with and without absorption respectively, (b) 
Position of poles and zeroes corresponding to curves A and C in the k 2 plane. When 
poles and zeroes lie very close to each other, the marker resembles an asterisk. 

Table 1. Poles and zeroes of the S matrix in k plane for the rectangular barrier of 

eq. (3). 

U B = B fm- 2 , a = 5 fm, b = 5.5 fm 



Poles 



Zeroes 



W = 0.2 fm 



-2 



W = 0.2 



0.582-2.100 x 10~ 3 J 
1.642-8.461 x 10- 3 i 
1.747-1.922 x 10- 2 i 
2.330-3.452 x 10' 2 / 
2.915-5.444 x 10~ 2 z 
3.500-7.908 x 10~ 2 i 



0.605-0. 167 i 
1.67 1-9.332 xlO~ 2 i 
1.747-7.535 x 10~ 2 i 
2.329-7.614 x 10~ 2 / 
2.914-8.735 x 10~ 2 / 
3.499-0. 1 06i 



0.582 + 2. 100 x 10- 3 z 
1.642 + 8.461 x 10~ 3 i 
1.747+ 1.922x 10~ 2 / 
2.330 + 3.451 x 10~ 2 / 
2.915 + 5.444 x 10~ 2 i 
3.500 + 7.908 x 10" 2 i 



0.605-0. 163 i 
1.675-7.641 x 10~ 2 i 
1.748-3.691 x 10- 2 / 
2.331-7.1 19 x 10- 3 i 
2.916 + 2.151 x 10~ 2 z 
3.501+ 5.197 x 10~ 2 i 



bottom of the pocket, the S matrix pole and zero are located very close to each other and 
they gradually diverge as far as the imaginary part is concerned, as the real part of the 

nole.s and 7P,roe.s in thf, p.np.rpv nlanp. inrrpasp.s Thus nnp. ran <stafp that np.aVc in /r , 
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Figure 2. (a) Variation of a & and <r a b with energy for the Gaussian barrier of 
eq. (4) with U B = 10 far 2 , W = 1 fm~ 2 , a = 0.5 fm, a/ = 0.3 fm, R o = 4.25 fm 
and R = 4.25 fm. Curve A is the plot of cr a b and curves C and B correspond to cr s 
with and without absorption respectively. The arrow mark indicate the position of 
the first extremely narrow peak in curve B which have not been reproduced due to 
the finite step size of the energy intervals, (b) Same as figure Ib for the potential 
of eq. (4). 



indicate the presence of resonance states generated by the complex potential. The 
conclusion is, non-elastic or reaction channels provide a much clearer signature of 
resonance than the corresponding elastic channel. Figure la also indicates that the 
width of the cr s resonance peaks in curve B gradually increases as the real part of 
the energy becomes large as it should since Foe Im(ZiR). Physically speaking, this 
means that a resonance state with higher energy can decay more easily by tunneling 
across the barrier as compared to the states with lower energy. The conclusion from 
figures la and Ib is that peaks in cr ab arise from the zeroes of the S matrix but at the same 
time indicate resonance state positions because the real part of zeroes and poles for these 
cases are practically the same as indicated in figure Ib. The first peak in cr ab at 
k 2 ~ 0.6fm~ 2 is less prominent because in this case, the S matrix pole and zero having 
appreciable imaginary parts are very close to each other thus making \S(k) ~ 1|. This 
point will be elaborated later. 
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.Poles Zeroes 

W = 1 frrT 2 W = OfrrT 2 W = I far 2 



0.787-0.416 x 10~ 4 z 
1.561-1.128 x 10- 3 i 
2.303-6.975 x 10~ 3 / 
3.001-3.520 x 10~ 2 z 
3.668-0.1 13x 


0.946-0.525/ 
1.591-0.309i 
2.312-0.214/ 
3.001-0.191/ 
3.664-0.242* 


0.788 + 1. 416 x 10~ 4 / 
1. 561 + 1. 128 x 10- 3 f 
2.303 + 6.975 x 10~ 3 / 
3.001 + 3.520 x 10~ 2 ; 
3.668 + 0.1 13i 


0.946-0.525/ 
1.591-0.307i 
2.314-0.201/ 
3. 008-0. 121 i 
3.676-1.510 x 10~ 2 / 



2.2 Gaussian barrier with absorptive pocket 

Now, we examine the same problem for a smoother potential having a Gaussian barrier 
with an absorptive Woods-Saxon pocket 



J r -^L\ 






(4) 



with U B = 10 far 2 , W = 1 far 2 , R G = 4.25 fm, R = 4.25 fm, a = 0.5 fm, a t = 0.3 fm. 
In this example, the barrier is comparatively broader. We have solved the Schrodinger 
equation using the Runge-Kutta method and the results of the calculation are summarized 
in figures 2a, b and table 2. We see that this potential essentially reproduces results that 
corroborate those of the previous potential, given by (3) leading to similar conclusions. 
We note that the sharp resonance peaks in a s in curve B for W = gradually gets 
narrower (~ 10~ 4 fm~ 2 ) until they are too narrow to be observed due to the finite size of 
the energy steps used for the calculation. The arrow in figure 2a indicates the position of 
the sharpest resonance of curve B. But as in the previous example, the locations of these 
extremely sharp peaks match with the positions of the corresponding absorption peaks for 
the higher lying states. This point is explained in 3. Compared with the former case, this 
potential has a larger imaginary part and width and the absorption cross section shows 
less structure as compared with that in figure la. Numerically speaking, this is due to the 
fact that poles and zeroes of the S matrix for this complex potential are much closer to 
each other than the corresponding case depicted in figure Ib. In this case, barrier 
tunneling is reduced due to the larger width and higher barrier. As expected, in the present 
case the resonance peaks in cr s for W = are sharper than the previous case studied. 

2.3 Woods Saxon type barrier with absorptive pocket 

As the last example, we analyse the potential having a still larger barrier with a Woods- 
Saxon shape and an absorptive pocket of the same shape [11]. This potential is given by 

, r< R(l 

(5) 



The parameters used for this potential are U B ~ lOfm 2 , W 1 f m 2 , RQ =4.5fm, 
fli = 0.12 fm, 2 = 0.25 fm, a w = 0.05 fm, c\ = l.Sfm, c 2 = 2fm, cw = l.Sfm. The 
scattering and absorption cross sections as well as the positions of the poles and zeroes 
given in figures 3 a, b and table 3 at first glance, do not seem to be similar to our previous 
two potentials of (3) and (4). The large barrier width effectively suppresses tunneling 
though the barrier and the corresponding resonances in curve B with W have 
extremely narrow widths. These peaks in curve B have not been reproduced in the graph 
due to the finite step size of the energy interval but we have indicated them by arrow 
marks. We have verified that these fine peaks in cr s are reproduced with discernible width 
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Figure 3. (a) Variation of a s and <r a b with energy for the Woods-Saxon shaped 
barrier of eq. (5) with U B = 10 far 2 , W = I fm" 2 , a]=0.12fm, a 2 = 0.25fm, 
a\>/ = 0.05 fm, c\ = l.Sfm, c 2 = 2fm, GW = l.Sfm, 7?o = 4.5 fm. Curve A is the plot 
of cr a b and curves C and B correspond to <j s with and without absorption respectively. 
The arrow marks indicate the positions of the extremely narrow peaks in curve B 
which have not been reproduced due to the finite step size of the energy intervals. 
Below E = 10fm~ 2 , curves B and C practically merge into one. (b) Same as figure Ib 



C2 = 2.0 fm, cw = 1.5 fm 



Poles Zeroes 



W = OfrrT 2 W = 1 fnT 2 W = OfirT 2 W = 1 far 2 

0.978-3.244 x 10' n i 1.081-0.459i 0.978 + 3.286 x l(T n i 1.081-0.459* 

1.928-7.743 x KT 9 z 1.945-0.246i 1.928 + 7.743 x 10- 9 j 1.945-0.246/ 

2.808-2.274 x KT 5 / 2.815-0.146/ 2.808 + 2.274 x I0~ 5 i 2.815-0.146J 

3.275-4.774 x lO' 2 ! 3.267-5.976 x lQ~ 2 i 3.275 + 4.774 x W' 2 i 3.283 + 3.719 x KT 2 / 

3.518-9.981 x 10- 2 i 3.504-0.147/ 3.518 + 9.981 x lO' 2 / 3.545 + 5.013 x 1(T 2 / 

3.779-0.180i 3.779-0.229J 3.779 + 0.180i 3.771 +0.128J 



if the width of the barrier is decreased. More importantly, the display of the zeroes and 
poles of the S matrix does not show the symmetry in the E-plane which was present in the 
earlier two examples. This indicates that the symmetry in figures Ib and 2b should not be 
considered as the general feature of the zeroes and poles of the 5" matrix for complex 
potentials. However, one feature which may be noted is that as far as the deeper lying 
resonances are concerned, both the zeroes and poles practically coincide with each other. 
In such cases, peaks in cr ab also are not observed. The nature of variation of poles and 
zeroes in figure 3b shows the dramatic difference between below barrier absorptive 
pocket resonances and above barrier resonances. The above barrier resonances in this 
case are very clearly manifested in cr ab but not in a s . It may be mentioned here that the 
resonances in heavy ion scattering which have highly absorptive pocket followed by large 
Coulomb barrier are often associated with barrier region orbiting resonance states. 

3. Discussion and conclusions 

We further analyse our results as follows. In the case of real potential we have 
Re(fcp) = Re(&z) and in the examples studied, we find that the same relation is valid to a 
good degree of accuracy for complex absorptive potentials. Since the poles are observed 
to occur lower down in the &-plane as compared to the corresponding zeroes, the 
following relations also hold for complex potentials. 



> lm(k P ] and |I 
For the s-wave scattering case, we have 



Near a pole or zero, we can write the S matrix as 



~ \S(k] f = and -^ \S(k) \ 2 > at k = Re (k P ) 

Q/C Q/C 

and hence S(k) has a minimum at k = Re (&/>). This will give rise to a maximum in the 
absorption cross section at the same energy. However, the corresponding rise in the 
scattering cross section is suppressed by the imaginary part of the potential due to the fact 
that in this case, 

Im(* 2 



Im (k 

and for all cases, either Im^z) ~ Im (kp] or |Im(&2:) < |Im(fc/>)|. Hence |1 S(k}\ A 
becomes much smaller than the maximum possible value" 4. On the other hand, in the case 
of real potentials, |1 S(k}\ = 4 near a resonance. Thus the lack of mirror symmetry 
with respect to the real k axis and the nature of the poles and zeroes for complex 
potentials suppresses the resonance peaks in cr s which the corresponding real potential 
can manifest when proper resolution is taken. 

The width AE of the absorption peak of (1 - |5(A:)| 2 ) can be shown to be 
AE 1 ~ 4\R.e(kp)lm(kp)\. Even though a similar relation is valid for the peak in cr s , 
due to the drastic suppression of the peak height as explained above, these are suppressed 
in <7 S . 

As seen earlier in the case of the potential given by eq. (5), there are no peaks in u ab 
observed at lower energies because in this region, the complex zeroes and poles having 
significant imaginary parts are practically same leading to \S(k}\ ~ 1 for E = E R . This 
can also be interpreted in terms of the inhibition posed by the large barrier to generate 
absorption in the pocket at lower energies. Thus in general, the variation of cr a b provides a 
better indicator of the effective resonances in complex potentials as compared to cr s . 

To summarize, we have examined the poles and zeroes of the S matrix for three typical 
complex potentials having absorptive pocket followed by a barrier. We find that in 
general, the real parts of the low lying zeroes and the corresponding poles of S(k) are 
practically the same and the imaginary parts are significantly different either in sign or 
magnitude or both. Such a statement is also substantially true for higher poles and zeroes. 
Due to this, |1 S(k}\ 2 gets greatly suppressed as compared with its peak value (= 4) 
and hence elastic channel resonance peaks are not manifested appreciably. On the other 
hand, the zeroes of the S matrix corresponding to a pole enhance the absorption cross 
section factor (1 |5(A;)| 2 ) in the vicinity of resonance energies. This enhancement 
prominently occurs when | Im (kp)\ is significantly larger than | Im (kz}\ but when Im (kp] 
is quite close to lm(kz], even in <r a b the peaks get eliminated. As a result, peaks in a & ^ can 
be correlated with the resonances generated by the complex potential. This is consistent 
with the fact that in the analysis of resonances in nucleus-nucleus collisions, one of the 
important criteria for identifying resonance is that resonance structure should appear 
correlated in energy in several exit i.e., reaction channels. [12] The peaks in the absorption 
cross section also gets suppressed if the width of the barrier is quite large and in such 
cases, peaks in cr a b occur due to the barrier top states. These observations mean that in 



be comparatively weaker. 
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1. Introduction 

A new irreversible thermodynamic formalism, viz. the generalized phenomenological 
irreversible thermodynamic theory (GPITT) [1-6] has been developed. This theory reveals 
that every irreversible thermodynamic formalism is necessarily governed by the validity 
of Gibbs-Duhem equation (GDE). Hence, according to GPITT all nonequilibrium thermo- 
dynamic situations belong to GDE-regulated regime. The breakdown of GDE means the 
crumbling down of irreversible thermodynamics itself. For example, GPITT derives (does 
not assume) the traditional local level version of the Gibbs relation, namely: 

dt dt dt ^ dt 

under the condition that its corresponding GDE, namely 
dr dp 



density, ^ is me cnemicai potential per unit mass or me component /e, Ck is me mass 
fraction of the component k and t is time. In GPITT one also derives an extended Gibbs 
relation, 



which is valid provided its GDE, namely: 



also holds true, where QJ and /3j are the extensive and intensive parameters respectively 
and are identified in GPITT as the additional irreversibility parameters. They originate 
from the additional nonequilibrium interactions [3-5] other than the thermal, mechanical 
(both long and short range ones) and chemical. Unfortunately, at present no physical basis 
for these so-called additional nonequilibriura interactions is in sight. 

Thus the classical irreversible thermodynamics (CIT) of Onsager-Prigogine-Meixner- 
de Groot [7-11] is reproduced in GPITT (c.f. eq. (1.1)) but with a difference that 
it stresses the fundamental role of the GDE, eq. (1.2) and one is not at all required 
to invoke the local equilibrium assumption for the derivation of eq. (1.1). On the 
other hand the derived (not proposed) extended Gibbs relation, eq. (1.3), of GPITT has 
specific requirements and those do not seem to be fulfilled by the existing other 
irreversible thermodynamic theories, for example extended irreversible thermodynamics 
(BIT) [12-15]. 

The strength of the GPITT lies in its simple but overall picture about the thermo- 
dynamics of irreversible processes, though at times a nonequilibrium situation may be a 
highly complex one from many points of views. Amongst the vital tenets of GPITT is the 
recognition of an universal inaccessibility principle (UIP), which reads as: 

"In every neighbourhood of every equilibrium and nonequilibrium state of a system 
there exist states (both equilibrium and nonequilibrium ones) which remain 
unattainable by reversible or spontaneous paths."" 

This principle guarantees the existence of an integrating denominator to the concept of 
thermodynamic heat also identified and defined so in GPITT. This then leads us to the 
corresponding entropy function and the thermodynamic temperature for nonequilibrium 
states. Indeed, while presenting our development of GPITT [1] we had simply stated the 
UIP and used it without elaborating its essential contents. However, recently we have 
demonstrated [3] that UIP serves as a local level analogue of Clausius' inequality. Thus, 
we take this opportunity to elaborate on its essential contents. 

The next section describes reversibility and spontaneity. In 3 some relevant concepts 
of GPITT are presented, in 4 and 5, how the use of UIP is involved in the De 
Donderian thermodynamic description of chemical reactions at finite rates and in the CIT, 
is demonstrated. Section 6 contains remarks on the inclusion of the situations with 
negative temnerafurfis in the fold of TJTP and 7 contains concluding remarks. 



In thermodynamics of equilibrium states, a prescribed path through the succession of 
equilibrium states is said to be reversible because the same succession of equilibrium 
states can be traversed in the reverse order. These paths in principle can only be traversed 
infinitely slowly and hence they are also termed as quasi-static paths. Recently, the 
reversible paths through the succession of nonequilibrium stationary states (NSS) have 
also been recognized. Keizer has developed a generalized fluctuation-dissipation 
principle, a statistical kinetic theory description, and one of its results is the prescription 
of reversible paths through NSS [16-20], which too by prescription can only be traversed 
infinitely slowly, meeting the requirement of quasi-staticity. Therefore, in this 
presentation reversibility implies the passages through the succession of either 
equilibrium states or NSS. However, along every reversible path through a succession 
of equilibrium states a set of corresponding time independent differential conditions, 
G/(y,-)s are followed, namely: 

' (2.1) 



where ;y,-s are the corresponding variables. Whereas for the reversible paths through NSS 
there exist similar equations, 

/(r) =J%,W)))- (2-2) 

To distinguish between equilibrium and nonequilibrium we have used x,s in eq. (2.2) 
instead of y f s and changed to a local level description. However, for the present purpose it 
is not required to know the exact expressions for the above functions. It is only sufficient 
to know that there exist the constitutive equations, gj(xi(r)}s and G/(;y,-)s. 

In the true sense of the word 'spontaneous' there are no processes which occur without 
intervention. But still one describes the so called spontaneous processes. These are the 
processes which cannot be arrested if the constraints on the system are not appropriately 
altered. The constraints on the system may be time invariant or time dependent. In UIP 
and GPITT the evolution of a system under the given constraints is described as 
spontaneous one. Let us make it a specific mathematical idea. Consider a closed system 
in an equilibrium. Now if the system is instantly brought into a nonequilibrium condition 
by a suitable manipulation it will follow an irreversible adiabatic path if the adiabaticity is 
adhered throughout. Depending upon the initial time invariant conditions to which the 
closed adiabatic system is subjected to, a final equilibrium state is reached. To reach 
different final equilibrium states adiabatically, the initial conditions need to be changed 
accordingly. This behaviour is well documented and is easy to visualize [21, 22]. Thus for 
different irreversible adiabatic transitions, even if all of them originate from the same 
initial equilibrium state by suitable manipulations, the end equilibrium states are different. 
That is, along a given irreversible adiabatic transition of a system it passes through an 
unique succession of local nonequilibrium conditions (within the system) or in other 
words with a given set of dynamical conditions there is one and only one real trajectory 
for an irreversible adiabatic transition. 

However, to ascertain which trajectory is followed during an irreversible adiabatic 



where g/s are the dynamical differential conditions or constitutive equations whose 
variables are the x,-(r,f)s and r is the position vector. Thus the knowledge of all the 
g/(.Xi(r, f))s along with the initial conditions identifies an irreversible adiabatic. However, 
every spontaneous transition whether it is an adiabatic or a nonadiabatic one possesses the 
same characteristics, that is, the passage is determined by gj(x((r, f))s. Therefore, eq. 
(2.3) contains a set of constitutive equations of spontaneous paths. The existence of the 
constitutive equations g/s in terms of jc/(r, t)s amounts to the recognition of the transitory 
local conditions as local nonequilibrium states. If a process is much faster than the time 
required to experimentally register z/(r,?)s then no g/s can be prescribed. Some of the 
examples of jt/(r, r) are temperature, pressure, concentrations, internal energy, volume, 
physical fluxes, gradients of intensive quantities etc. Recall that even in nonequilibrium 
situations one measures at the local level temperature, pressure (or stresses) and 
concentrations with ease and also the physical fluxes. For example, temperature can be 
measured within the time period approaching nanoseconds and in situ measurements of 
concentrations require a time period much shorter than a picosecond (for example the 
spectrophotometric absorption and emission measurements). Thus with the recognition of 
the constitutive equations a notion of nonequilibrium states gets informally incorporated 
in our description and irreversible thermodynamics makes it a formal concept. From the 
above discussion it is clear that by definition the local nonequilibrium states exist for 
exceedingly short but finite duration depending upon the fastness of the process. Further, 
it is not necessary that all the Xf(t)s would become independent thermodynamic variables. 
The latter are determined by the laws of thermodynamics and not by a constitutive theory 
[3, 5]. In nutshell g/(;c/(r, r))s and jt/(r,r) = jc,-(r, 0) determine which spontaneous path 
is being followed by the system. These we know say from fluid dynamics. However, 
in connection with the application of UIP we will see in the following that all g/s are 
not required to construct the corresponding condition, namely /(r, t). Also notice that 
eqs (2.1) to (2.3) are the constitutive theory. 

To cite the illustrative examples of the above deductions in terms of processes we 
may recall that under a given condition a chemical reaction proceeds through a 
corresponding given succession of nonequilibrium states. If this were not true then it 
could not have been possible to obtain the reproducible values of rate constants on 
repeating the experiment. Recall from the literature on chemical reaction rates that 
even in the case of diffusion controlled reactions the experimentally measured accurate 
values of rate constants are available. Yet another simple example is that of the rigid 
body heat conduction. When a heat conducting bar is placed between two heat reservoirs 
at different temperatures and if the temperature at a position z along the length of 
the bar is measured with time then the same temperature-time profile is obtained on 
repeating the experiment. A good example of this is the experimental recording of 
temperature waves [23]. To this list one can add numerous examples for there is no dearth 
of them. 

Further, recall that in the case of irreversibility the end equilibrium states are mutually 
inaccessible; the accessibility is in one direction only. For example, from an equilibrium 



local nonequilibrium states cannot be traversed in the reverse direction. In contrast 
Muschik [21] describes that any two equilibrium states of a closed system are attainable 
by more than one irreversible adiabatic paths but then it demands at least one additional 
equilibrium state to be reached in between [21]. Therefore, the contents of UIP should not 
be confused with such situations. Yet another example, if sufficient care is taken, it 
can mislead to that of a reversible chemical reaction. Such reactions can be made to 
proceed in forward as well as in reverse directions so that the chemical composition of the 
system is changed exactly in the mutually reverse order in the two situations. But if 
the chemical conversion is not carried under the reversible conditions then one needs 
to adequately change say T and/or p to effect the chemical reaction at a finite rate 
in the reverse direction and hence although the chemical composition changes exactly in 
reverse order it is under a different thermodynamic condition, that is the system traverses 
a different succession of nonequilibrium states. Thus the above elaboration makes it 
precise what is implied by the terms spontaneity and reversibility in UIP and hence in 
GPITT. 

Moreover, all the heat differentials used in global and local level descriptions, in the 
following sections, are the inexact differentials but they have not been distinguished so in 
their notation compared to the other thermodynamic exact differentials. Hence the former 
should not be read as a differentiation of a function. Further, in order to appreciate UIP 
and its application we felt that it is advantageous to describe first the appropriate concepts 
of GPITT in adequate details which is the subject matter of the next section. 



3. Some relevant concepts of GPITT and application of UIP 

To meet the requirements of self-sufficiency and brevity we are describing below the 
thermodynamic heat as defined in GPITT [1-5] considering only a simplest situation. Let 
the following fluid-dynamical internal energy balance equation (which does have suffi- 
cient element of generality) be a total description of a situation in hand: 

p~= -divq + T: gradu. (3.1) 

at 

In eq. (3.1) q is the so-called heat flux density, T is the stress tensor and u is the bary- 
centric velocity. 

In the first step the stress tensor is split into two parts, 

T = -/?l + n, (3.2) 

such that the first term on the right hand side of the following equation, 

dv 

T: gradu = pp + II: gradu (3-3) 

at 

completely quantifies the nondissipative part of the work due to the contact forces. 
Therefore, the second term on the right hand side of eq. (3.3) measures the total 



are designated as Born s neat, 

p ~- = -div q + FT. grad u. (3 .4) 

dt 

Note that according to GPITT eq. (3.4) is the definition of Bern's heat corresponding to 
the internal energy balance equation of eq. (3.1). 

Since even in CIT [7-1 1] the terms on the right hand side of eq. (3.4) taken together are 
not adequate enough to describe the total change in entropy of the system we must search 
for a new heat concept. For a general nonequilibrium situation GPITT defines the 
thermodynamic heat as 

dq dq E dq nc 



Note that the nonconventional (nc) heat is not the one designated in the literature as 
Clausius' uncompensated heat, though through the former one can arrive at the corres- 
ponding expression for the latter [3]. Notice that the nonconventional heat is also recog- 
nized for the first time in GPITT [1-5]. In support of the above statement it can be shown 
that in the thermodynamics of equilibrium states of open systems the equilibrium Gibbs 
relation follows elegantly by using the nonconventional heat [24]. 

With the combination of eqs (3.1) and (3.4) and the standard definition of enthalpy, h 
(per unit mass), 

h = u +pv (3-6) 

one gets the following expression for Bern's heat: 

^B_^ ,-i d ^ (37^ 

dt ~dt~ P (3>7) 

which resembles with the corresponding equation in the thermodynamics of equilibrium 
states and hence the first law of thermodynamics in the local form, eq. (3.1), is obeyed by 
Bern's heat only. 

Now we are in a position to use the UIP. From the discussion of 2 above, it follows 
that with the restriction of spontaneity, that is by a given equation of a spontaneous path, 

^(r,r)=/(r,r) (3.8) 

there remain numerous nonequilibrium and equilibrium states of the system unattainable. 
This is what UIP precisely describes. Notice that only relevant dynamical conditions 
gj(x{(r, r))s need to be incorporated to construct /(r, t] of eq. (3.8). Therefore, from the 
inaccessibility described on following (3.8) it immediately follows that dq has an 
integrating denominator and thus we obtain the entropy (s) and the local thermodynamic 
temperature (T) in a nonequilibrium situation as 

ds - l dq (3 9^ 

-' (3 ' 9) 



the name entropy for the function s in contrast to the practice in other irreversible 
thermodynamical theories where one encounters the terms like generalized entropy, 
calortropy and an entropy like function [12-14]. The remaining parts of the details 
wherein the traditional Gibbs relation, eq. (1.1) and the extended Gibbs relation of 
GPITT, eq. (1.3) are arrived at are to be found in our earlier discussions [1-5]. 

On the other hand, on following the restriction of reversibility (for closed systems), that 
is the condition of eq. (2.1), 

<1Q C = F, (3.10) 

where dQ c is the reversible differential element of heat and the subscript c stresses the 
conduction mode of energy exchange, no nonequilibrium state is attainable from any one 
of the equilibrium states of the system. Thus we conclude that all equilibrium states of a 
system are mutually accessible. This is the equilibrium part of UIP. In this case the form 
of the function F is none else than the first law of thermodynamics which in the simplest 
case reads as F = dU + pdV. Moreover, on using F = 0, the differential condition, in eq. 
(3.10) (that is on imposing the condition of adiabaticity) the Caratheodory principle is 
obtained which describes that not only no nonequilibrium states are attainable but 
numerous other equilibrium states too remain unattainable on following a reversible 
adiabatic path. Notice that traditionally while applying Caratheodory 's principle no 
nonequilibrium states are under consideration the unattainability of the latter category of 
states is neither of any consequence nor is required in thermodynamics of equilibrium 
states. The case of reversible paths through NSS is considered in one of the following 
sections. In order to make it further clear that UIP is not an offbeat statement the next two 
sections demonstrate the implied involvement of UIP in the already existing irreversible 
thermodynamic descriptions. For this purpose some relevant stepwise details of the two 
chosen irreversible thermodynamic descriptions are described in the next two sections in 
such a way that a connection with the various concepts of heat of GPITT is established. 

4. The derivation of a variance of UIP from the De Donderian thermodynamic 
description of chemical reactions at finite rates 

Though the Belgian group [7,24] does not seem to claim it emphatically, the De 
Donderian thermodynamic description of chemical reactions within a closed system does 
constitute one of the first irreversible thermodynamic descriptions. Prigogine and Defay 
in their classic monograph on chemical thermodynamics [25] describe it very elegantly. 
They start with a closed system having irreversibility within it only due to the occurence 
of a chemical reaction. Since the system considered is a closed one, the following form of 
Clausius' inequality is used: 

as>, (4.1) 



where dQ c appears in the corresponding first law expression, namely: 



energy and V is the volume of the system. Since no irreversibility is assumed in heat and 
volume exchanges the temperature 7\ and pressure p, are those of the system itself. In the 
next step eq. (4.1) is converted into an equality by identifying a quantity of heat dQ', 
termed as the Clausius uncompensated heat having the property, 

d<2' > 0, (4.3) 

which vanishes in equilibrium and hence in the limit of reversibility through the 
succession of equilibrium states. Indeed, barring a sole example [26], hardly one finds the 
latter assertion in the literature. In view of eq. (4.3) the Clausius inequality, eq. (4.1), 
produces the following equation for the entropy change: 

d, = ^M. (4.4) 



Since, the uncompensated heat measures the irreversibility on account of the occurence of 
a chemical reaction, it is then quantified as 

d<2' = Ad > 0, (4.5) 

where A is identified as the chemical affinity and is the extent of chemical reaction. 
Next eqs (4.2), (4.4) and (4.5) are combined to give 



dU = TdS - pdV - Ad. (4.6) 

We notice that an irreversibility due to a chemical reaction can arise only when it 
occurs at a finite rate, this is what the Belgian school failed to stress. Since by definition 
A vanishes at equilibrium it needs to behave in identically the same way along a reversible 
path even if during a reversible transition d 7^ 0. From which it immediately follows in 
the view of the definitions of eqs (4.4) and (4.5) that along a reversible adiabatic of a 
closed system even if chemical composition within it changes on account of accompanied 
internal chemical conversion the reversible adiabatic remains perfectly isentropic. We 
have demonstrated this in one of the previous discussions by using the flux-force 
relationships for chemical reactions [5]. It is, therefore, surprising to learn that this simple 
looking deduction failed to catch an attention in thermodynamic literature so far. 

We stress that eq. (4.6) is, indeed in the Pfaffian exact differential form but it is at the 
same time a nonequilibrium thermodynamic description too. In other words, the above 
descriptions should involve the time parameter, so that eq. (4.6) reads as 

W T f-p%-A* (4.7) 

dt dt dt dt v ' 

That is we have 

U=U(t) t S=S(t), V=V(t), f=f(0, T=T(t), p=p(0, A=A(t}. (4.8) 
For- the same reason eq. (4.4) should read as 

dS 1 \dQ c dQ'] 

= - + . (49) 

At T\ At ^ At \ ^ y ) 



Clausius uncompensated heat, because the calorific quantity dQ c of eq. (4.2) is inade- 
quate to define the total change of entropy. Of course, dQ' quantifies the amount of non- 
retrievable damage caused by a chemical reaction at a finite rate. 

Now let us view the above equations in terms of the concepts contained in GPITT. For 
this purpose we denote by dQ, the sum 

dQ = dQ c + dQ' (4.10) 

and call it the thermodynamic heat. The general definition of thermodynamic heat is 
given in GPITT [1-5] as in eq. (3.5) above. The name thermodynamic heat stems from the 
fact that it directly leads us to thermodynamic thermal parameters S and T. Thus eq. (4.9) 
is compactly expressed as 

dS ld <2 



The acceptance of the definition contained in eq. (4.1 1) implies that in this nonequilibrium 
situation the thermodynamic temperature has been identified in T and it has been assigned 
as an integrating denominator to the heat differential dQ. Indeed, T so obtained is a 
temperature corresponding to the nonequilibrium state of the entire chemically reacting 
system, because the system is assumed to be thermally homogeneous. Hence, 5 so 
obtained, is the entropy in a nonequilibrium state originating due to a chemical reaction at 
a finite rate. 

From Caratheodory's theorem [22, 27] we know that the existence of an integrating 
factor for an inexact differential implies the existence of a differential equation, in terms 
of the same inexact differential, describing the corresponding unattainability. In the view 
of this and eq. (4.11) the following statement immediately follows: 

"dQ is an inexact differential but is integrable, that is it possesses an integrating 
factor, which implies that there exist states which remain unattainable by the paths 
given by dQ/dt = F(t)." 

The function F(t) is defined on the lines of eq. (2.3) and on keeping in mind that there is 
no spatial inhomogenieties within the system under consideration. Thus on using the 
above condition, 

^ = F(t} (4.12) 

in eq. (4.11) we have 

dS = ^F(t}dt (4.13) 

whose solution (in view of eq. (4.7) and in general T will be time dependent) reads as 
S = S(Z/(f),V(0,(0) (4-14) 

subject to an additive constant. 



= Jr(t) 

where dQ is the thermodynamic heat defined in equation (4.10). 

Thus UTP-I is nothing else than the UIP applied to a global nonequilibrium system 
having thermal and mechanical homogeneity. Notice that in this case we have 
F(t) = dU/dt + pdV/dt + dQ'/dt, Moreover, on making the dynamical condition F(t] 
of eq. (4.12) as a time independent constant, the entropy varies linearly with time and on 
using F(t] = in eq. (4.12) the system gets relegated to the corresponding NSS. But both 
the cases cannot be realized in a real situation for a chemically reacting closed system. 
Hence the case of F(t) = means the attainment of chemical equilibrium and F(t) equal 
to a time independent constant is devoid of a physical significance in the present case. 

Note that if there exist irreversibilities for example due to the thermal and mechanical 
inhomogeneities within the system the above global inaccessibility principle, UIP-I, 
cannot be applied as such. For such situations a local level consideration is demanded. 
This we describe in the following section. 

5. UIP from CIT 

In this and in the following sections all the per unit mass time derivatives are in fact the 
substantial time derivatives at a given value of the position vector, r. We assume that the 
illustrative system at hand is adequately described by the internal energy balance eq. (3.1). 
Note that in eq. (3.1) we have w(r, t), p(r, f), q(r, f), T(r, t) and u(r, t). However, if the r 
and t dependence has not been depicted for the quantities in the following expressions of 
this paper it is implied as usual. In the. development adopted in CIT [7-1 1] the stress 
tensor is splitted as 

T = -/>1-P, (5.1) 

where 1 is the unit tensor, P is the viscous pressure tensor and p is understood as the 
hydrostatic pressure or more precisely -(l/3)TrT. Next the validity of the traditional 
Gibbs relation, eq. (1.1), is taken as a thermodynamic condition of the existence of the 
local equilibrium within the system. Then (3.1), (5.1) and (1.1) are combined to give the 
following expression for the rate of change of the local entropy, 




dt 

\ K / j " r K 

1 A HA 

(5.2) 
In arriving at eq. (5.2) the following standard fluid-dynamical equation is used: 

= -div T + m (53} 

dt ' 

where 3k are the diffusion fluxes, m^ are the rates of mass fraction changes on account of 

rhpmirfi] rpar-tinn \c fViA r-nrrpcnrnHincr pYt^nf rf rhpmirfll reaotinn anH tine*. 



where M* is the mass per mole of the component k and i/*s are the stoichiometric 
numbers. 

Now we proceed to establish a connection between the above CIT description with the 
concepts contained in GPITT as follows. Recall from GPITT the definition of Bom's 
heat, dq%, which in the present case is quantified as 

p : vu (5.5) 



. 

at 

and recast eq. (5.2) as 

pT -/= -divq-P:Vu + y /4fcdivJ*4--4-r- (5.6) 

at *-* at 

Thus the first two terms on the right hand side of eq. (5.6) constitute Bern's heat (c.f. 
eq. (5.5)) and are provided by the internal energy balance equation, eq. (3.1) hence they 
are calorific in nature. On the other hand, the last two terms of eq. (5.6) are the additional 
ones which are required for defining entropy rate and in GPITT are termed as 
nonconventional (nc) heat terms [3]. Therefore, in the present case we have: 



* Qt 

Thus in terms of Bern's and nonconventional heats eq. (5.2) gives 
ds 



dt T[dt dt 
The thermodynamic heat, dq, in this case (c.f. eq. (3.5)) is given by 

dq = dq E + dq nc . (5.9) 

Therefore, eq. (5.8) reads as 

* = I^. (5.10) 

dt T dt ^ ' 

Thus we see that the adoption of the traditional Gibbs relation, eq. (1.1) and the defini- 
tions of dqs and dq nc tantamounts to prescribing the definition contained in eq. (5.10) for 
the entropy rate. Therefore, it immediately follows that T should be identified as the local 
thermodynamic temperature in a nonequilibrium state and taken as the integrating 
denominator to dq, the thermodynamic heat element. Since the validity of eq. (5.2) is 
beyond doubt, the adoption of eq. (5.10) is implied therein. 

The assertions from eq. (5.10) are similar to those that follow after eq. (4.11) in the 
preceding section. The first statement from the above delineation is obtained as 

"d# is an inexact differential but is integrable, that is it possesses an integrating 
factor, which implies that there exist states which remain unattainable by the paths 
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,, 
in combination with eq. (5.10) gives 



(5-12) 



In this case, the local level inaccessibility principle is obtained from the above 
considerations as 

UIP-II: "In every neighbourhood of every local nonequilibrium state of a system 
there exist local states which remain unattainable on following the 
equation (dq/dt}(r,t] =/(r, r)" 

where dq is the thermodynamic heat defined in eq. (5.9). 

The discussion of 4 and 5 clearly surface out the operation of the corresponding 
inaccessibility principles viz. UIP-I and UIP-II. This has never been realized before. 
Recall that dQ/dt F(t] and dq/dt=f(r,t) are the mathematical statements of 
spontaneity as defined earlier in 2 and hence UIP-I and UTP-n are nothing else than the 
LOP stated in 1 above. We see that the form of /(r, t) in eq. (5.11) is/(r, r) = dw/dH- 
pdv/dt + dq nc /dt. Further, on using the restriction /(r, t) = in eq. (5.11) we relegate 
the system to NSS and thus UIP generates the reversible paths through the succession of 
NSS and thereby confirms the earlier proposal of Keizer [16-20] on a generalized 
thermodynamic footing. The use of /(r, t) equal to a time independent constant in 
eq. (5.11) means entropy varies linearly with time. Such situations could be realized in 
practice by appropriate manipulations of a system. However, it would be curious to know 
whether such situations are encountered or not in industrial processes and in biological 
systems. 



6. UIP and negative temperatures 

The general statement of UIP reproduced in the introductory section of this paper is to be 
understood only in terms of reversibility and spontaneity described in 2 of this paper. 
Thus the inaccessible states, under the restrictions of reversibility, are the nonequilibrium 
ones (from equilibrium states). While for the nonequilibrium states there are numerous 
nonequilibrium and equilibrium states in their vicinity which remain inaccessible by a 
given restriction describing the spontaneous path. Basically the nonequilibrium part of 
the UIP is indeed, a local level assertion and from which the corresponding global level 
assertions follow. 

In the preceding sections we have seen that UIP is adequate enough to arrive at the 
thermodynamic local temperature and the corresponding entropy function for a system in 
nonequilibrium. It is now our choice to coincide the temperature given by UIP with the 
Kelvin temperature and in this way we make entropy S or s so obtained a physical 
quantity exactly on the same lines as one does in the case of equilibrium. The physical 
realirv of the fhermometfic temnerahire for noneouiHbrinm states has been desr,rihp.rl 



temperature is the identically same physical entity irrespective of whether the system is in 
an equilibrium or in a nonequilibrium state. 

The negative temperatures are associated to the population inverted states and hence 
these states basically constitute the extreme levels of nonequilibrium. Note that, since no 
restrictions were placed on the system while defining the thermodynamic heat function 
and neither UIP involves any such restriction, hence the nonequilibrium states with 
negative temperatures also find themselves into the fold of UIP and its application. Note 
that UIP provides the function T and not a particular sign to it. Hence, it is for us to 
conform the T function so obtained, that is from UIP, to the conventions we follow for the 
temperature function. Thus the concern spelled out by Muschik [28] that the traditional 
inaccessibility principle of Caratheodory does not cover the states with negative 
temperatures is non-existent with UIP. 



7. Concluding remarks 

This presentation was aimed at the elucidation of UIP and to surface out the hidden use of 
it, which exists in the appropriate versions in various irreversible thermodynamic 
descriptions. It was felt advantageous to present this discussion in the manner of a gradual 
development of UIP and its application tracing the same right from the earliest 
nonequilibrium thermodynamic descriptions up to the present level. However, we have 
restricted ourselves only to the De Donderian thermodynamic description of chemical 
reactions at finite rates and CIT and we found that it is adequate enough to reveal the 
basis and the essential contents of UIP. In doing so we have proposed that the UIP 
based local thermodynamic temperature at ease can be coincided with the experi- 
mentally measured local temperature, of a nonequilibrium system, on the Kelvin scale 
and thereby the corresponding entropy function is obtained with the same footing as the 
equilibrium based entropy function albeit the former pertains to a nonequilibrium state. In 
passing we deduced that the reversible adiabatic transitions of chemically reactive closed 
systems remain perfectly isentropic even if the chemical composition within the system 
changes on account of concomitant internal chemical conversion. Further, it also gets 
revealed that the UIP does cover in its fold the nonequilibrium states of negative 
temperatures. Moreover, through UIP we also arrive at the reversible paths through the 
succession of NSS and thus confirm the earlier proposal of Keizer [16-20] in a most 
general setting. 

We have shown above that the Caratheodory principle is indeed a particular case of 
UIP which demonstrates the self-consistency of the latter. However, details of the 
correspondence between the two and related discussion is the subject matter of one of our 
forthcoming communications. Of course, Caratheodory's principle is the most economic 
statement of the second law of thermodynamics and in this respect nothing can surpass it 
but then it is also true that it is not capable of revealing all the details of the 
nonequilibrium behaviour of thermodynamic systems nor the entropy function for a 
nonequilibrium state can be generated through it. Thus the sacrificing of the said 
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1. Introduction 

The recent decades have witnessed the evolution of extended irreversible thermo- 
dynamics (HIT) [1-7]. It is claimed by the proponents of HIT that the classical irreversible 
thermodynamics (CIT) of Onsager-Prigogine-Meixner-de Groot [8-13] is not adequate 
enough to encompass all nonequilibrium thermodynamic situations. From the mono- 
graphs and review articles on HIT [1-7] one gathers a long list of phenomena suggested 
as not falling within the fold of CIT. A preliminary as well as a striking example quoted is 
that of temperature waves. A simple mathematical description of the wave propagation of 
heat is the Maxwell-Cattaneo-Vemotte (MCV) equation [14], namely 

q = -r^-AVr, (1) 

where q is the heat flux density, r is the relaxation time, A is the coefficient of thermal 
conduction, T is the temperature and t is time. Notice that the MCV equation is a 
constitutive equation and not a thermodynamic one. However, in EIT a thermodynamic 
derivation of it is presented and then it is claimed as a justification of the used thermo- 
dynamic framework. This approach of EIT shakes its very thermodynamic basis. 



equuionum or nonequuiDnum siaies ana corresponding processes, maeeu, n is IIOL ai an a 
business of thermodynamics to derive a constitutive equation [15-17], it simply takes the 
latter as it is given and presents the corresponding thermodynamic deductions. For 
example, thermodynamics cannot derive the ideal gas and van der Waal equations, they 
are taken in thermodynamics as they read. Yet another question one needs to ask is can a 
given constitutive equation lead us in constructing a suitable thermodynamic framework? 
The answer is yes and no. Yes, if the corresponding constitutive theory unambiguously 
establishes some additional fundamental interactions otherwise the answer is a clear no. 
In spite of this, for example for the phenomena of heat conduction BIT projects a picture 
that the MCV equation demands an extended thermodynamic framework whereas the 
Fourier law, namely 



q = -AVr (2) 

is consistent with CIT. The latter is also referred to as linear irreversible thermodynamics 
(LIT) because it is believed that the linear relationships between fluxes and forces such as 
eq. (2) are implied therein. This projection of HIT is an unfortunate situation in which 
thermodynamics remains guided by the individual constitutive equations and hence 
thermodynamics looses its fragrance of generality. 

Therefore, in this presentation we clarify this fundamental thermodynamic aspect by 
first considering some of the relevant features of the rigid body heat conduction followed 
by a generalized discussion using the framework of generalized phenornenological 
irreversible thermodynamic theory (GPITT) [18-22] which is a recent addition to the 
thermodynamic literature. The reason for the choice of the GPITT framework lies in the 
fact that in GPITT an extended Gibbs relation is derived (not proposed) using the laws of 
thermodynamics and hence it has specific thermodynamic requirements [20,21]. In view 
of this it is now expedient to analyse if the physical fluxes such as heat flux density etc. 
are taken as additional independent thermodynamic variables, as is the practice in HIT, 
whether the resultant thermodynamic structure would meet the thermodynamic 
requirements or not. Earlier Garcia-Colin et al [23] on the basis of kinetic theory 
concluded that the physical fluxes are irrelevant for constructing an extended thermo- 
dynamic framework. This presentation also deduce the same albeit on a generalized 
footing. Then, we proceed to describe some of the existing inconsistencies in most of the 
versions of HIT in order to shed some light on the existing mix-ups in them. Finally, we 
clarify that the experimentally measured temperature of a system in nonequilibrium is 
not a different physical entity than that of an equilibrium system. This presentation at 
its certain parts takes up some elementary level discussions too but for the sake of 
maintaining a clarity of thoughts. 



2. Some relevant features of rigid body heat conduction 

Consider a single component heat conducting rigid bar of uniform cross-section placed 
between two heat reservoirs of temperature T\(t] and T r (t). The subscripts 1 and r stands 
for left and right respectively. Let us assume that T\(t) > T r (t\. In this case we have the 



;re u is the internal energy and p is the mass density. Now T\(f) and T r (t) are mani- 
ated in such a way that at every position x of the bar the divq is time independent that 
we have, 

div q = const. (4) 

ice that q will be position and time dependent and const, in general will be position 
endent along the length of the bar. However, in view of (4), eq. (3) gives the following 
Lily of solutions, 

AM = const. Af . (5) 

; value of the const, of eq. (5) at a given position will of course depend also on the 
ditions imposed by the heat reservoirs. 

: rom eq. (5) we see that at a given position the same lapses of time bring about the 
ic amount of change in internal energy irrespective of the change in the value of q at 
: position with time. Thus the existence of a heat flux does not influence the internal 
rgy but its divergence does. It is understandable because the divergence quantifies the 
rejection or retention of energy. 

f const. = then the bar attains a steady state condition, q becomes position and time 
ependent and u becomes time independent but remains position dependent. Again the 
lie of u is obtained as position dependent in spite of the same value of q at all positions, 
his condition too the mere existence of q does not influence u. However, the stationary 
.es can be achieved for a series of different values of q by suitably manipulating T\ and 
Accordingly, at a given position the value of u will be different for different values of 
irough div q = in different set-ups. 

rhus on the one hand we see that q does not determine the value of u and on the other 
id by keeping say q constant u cannot be varied at a given position neither in a steady 
:e nor when div q ^ 0. Moreover, when q is allowed to relax to a value of zero or to a 
idy state value u does not remain constant. Thus even in this extremely fast process 
; relaxation time, r ~ 10~ 12 s [1]) u and q do not vary independent of each other. This 
icates that the mere presence of q cannot be taken as influencing the thermodynamic 
:e of a system. 

General considerations based on GPITT 

ITT constructs its framework by using zeroth, first and second laws of 
rmodynamics and hence it has a sound thermodynamic basis [20,21]. In contrast 
h BIT, GPITT is not aimed at the augmentation of LIT. GPITT yields a most general 
cription and from which LIT is obtained as a limiting case. One of the important 
elations of GPITT is that the Gibbs relation or an extended Gibbs relation becomes a 
rmodynamic description only by the virtue of the validity of the corresponding Gibbs- 
hem equations (GDE). That is, the crucial role of GDE is to confine our description 
11 within the thermodynamic domain. 
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~dt~T~dt + Tdt 4-T dr + r^' U dr ' (6) 

K. J 

where s is the per unit mass entropy, T is the local thermometric temperature, u is the per 
unit mass internal energy, p is local pressure, v is the specific volume, /J* is the chemical 
potential per unit mass of the component k, Q is the mass fraction of the component k, j3j 
and OLJ are the additional irreversibility parameters intensive and extensive respectively 
and O denotes the general scalar product. All the thermodynamic functions of me above 
equation are position (r) and time (?) dependent and for the sake of brevity we are not 
showing this dependence in the major part of this presentation. The derivation of eq. (6) 
involves the following standard fluid-dynamical equations, namely: 

p-^ = -divq-pp~ + II: Vu, (7) 

p =-div3k + pM k v k , (8) 

where q is the so-called heat flux density, 3k is the matter diffusion flux density of the 
component k, II is the dissipative stress tensor, u is the barycentric velocity vector, M k is 
the molar mass number, p is the mass density and is the extent of chemical reaction. 
When w, v and c k become time independent then according to eqs (7) and (8) q, II, Vu, 
3k, and p too become time independent and hence with respect to the three fundamental 
interactions, namely thermal, mechanical and those causing the changes in chemical 
composition, the system gets relegated to a stationary state. However, in this situation eq. 
(6) reads as 



Equation (9) describes that in the stationary state as determinable through eqs (7) and 
(8) the entropy of the system may still vary provided ajs((3js} remain time dependent. But 
the constancy of the above stated three fundamental interactions in a stationary state 
described through eqs (7) and (8) implies that q/ and /3y parameters must describe other 
than these three fundamental interactions which in GPITT are termed as additional 
nonequilibrium interactions [20-22]. This explains why in GPITT a/ and /3) parameters 
are termed as additional irreversibility parameters. However, at present nothing is known 
about the so-called additional nonequilibrium interactions. In view of the above neither 
the physical fluxes, namely q, II, 3k etc. nor the physical forces, namely VT, Vu, 
V(nk/T], etc. can be taken as a/ and/ or /?, variables. This assertion we further elaborate 
below by taking the case of rigid body heat conduction. 

Now again consider a single component heat conducting rigid bar of uniform cross- 
section. For the time being pretend that the heat conduction process can generate the so- 
called additional nonequilibrium interactions and for this purpose at least one pair of a 
and /? parameters would be required. In this case eq. (6) reduces to 

ds __ i d _ i _ da ... 



s s(u, a) = s(u(x, t),a(x, ?)) (11) 

and hence 

u = u(s,ai)=u(s(x t t),a(x,t)), (12) 

where x is the position coordinate along the length of the bar. We can as usual replace s by 
T in eq. (12) which gives 

u = u(T,a) = u(T(x,t),a(x,t)). (13) 

In view of eq. (13) the total differential of u reads as 
du dT du da 



where c v is the so-called heat capacity defined in this extended thermodynamic frame- 
work as 

- () 

According to eq. (3) in the case of a steady state heat conduction, that is q time 
independent and div q = 0, u too becomes time independent. However, in the case of 
steady state heat conduction our observations tell that T(x) also remain time independent 
and hence eq. (14) makes a(x) too a time independent parameter. Thus we conclude that 
T and a cannot be taken as mutually independent variables if these two parameters are 
considered as concerned with the phenomenon of heat conduction. Hence in the case of a 
rigid body heat conduction a and fi variables are redundant. The necessity of ex. and (3 
variables could be felt only when we also have other than the heat conduction process. 
This then clearly establishes that the heat flux q cannot serve as a variable in eqs (10) and 
(14). Therefore, for the case under consideration the appropriate equations, instead of 
eqs (10) and (14) are 

"- 1 - <> 



Now c v of eq. (17) has the standard definition, namely, 



instead of the definition of eq. (15). Thus in general eq. (17) rigorously gives 

u = u(T} = u(T(x,t}} (19) 

and hence from ea. (\f\\ we have. 



parameter, namely, the local temperature 1 , wnicn is none else man me local mermo- 
metric temperature. 

4. Some of the incoherences in EIT 

We point out herein two inconsistencies existing in EIT while it describes the phenomena 
of rigid body heat conduction. The first inconsistency is as follows. The EIT expression in 
the present case is [1-7] 

d77 = <T 1 dw-r 1 a-dq, (21) 

where 77 is said to be the nonequilibrium entropy, is the nonequilibrium thermodynamic 
temperature and is claimed to be different than the experimentally measured temperature 
T in the same nonequilibrium situation [1-4] and a is a parameter conjugate to the 
variable q. Equation (21) gives 

? = ?(,q), (22) 

that is 77 is considered as dependent on two independent variables. This then implies that u 
too becomes dependent on two independent variables say as u = u(rj, q) which may be 
expressed as u = u(T, q). Hence, the total differential of u should read as 

du _ du dT du dq 

cu = ar~d7 + aq'~d7' ( ' 

But in this connection in EIT eq. (3) is manipulated as [1-7] 

Jnn 



. (24) 

This means the function u and hence eq. (3) is not treated in the extended thermodynamic 
framework but for the 77 function eq. (21) is used. Thus what one is leading to is a mixture 
of unconnected mathematical descriptions which is not of a hybrid type though the latter 
if obtained cannot be termed as thermodynamic one. Strictly speaking one should use 
u u(Q, q) instead of u = u(T, q) but then eq. (24) would not be a direct result apart 
from the missing term corresponding to the variable q. However, this does not lead us any 
where. 

Secondly, there prevails another incoherence in the use of eq. (3) in EIT [1-7], It is 
rightly argued that u, the internal energy, can remain time independent provided div q is 
maintained effectively equal to zero. But it is not at all true that in general for a given set 
of different values of q there could be identically same constant value of u in the event of 
div q = 0. Equation pdu/dt = div q = gives a family of solutions of constant internal 
energy, namely, 

u = u(T] = u(T(x}} = const. (25) 

and not a single solution of constant internal energy. Let us elaborate it a little further. 

/-\ /vrt "f. y-k^-k*-i it* f^f-** * **! *-""Jl " "-C 4-* 1 A 



dx 
U= u(x)dx = u(x(T}} dT, (26) 



where u(x) is the internal energy per unit length of the bar at the position x. Now let us 
change Ti and T r to T{ and T( (T( > T' r }. In this new situation and under the corresponding 
condition of steady state heat conduction we have 

t L f T < dx 

U'= u'(x)dx = u '(x(T')}dT'. (27) 



o 

Obviously the temperature profile along the length of the bar will not be identically same 
in the two cases and hence U ^ U' implying u(x) ^ u'(x). Notice that in both the situa- 
tions we rigorously have div q = = div q at all the positions along the length of the bar 
but in general q' ^ q. This remains true even if the respective Ts are infinitesirnally away 
from T's. In this latter limit eqs (26) and (27) gives 



rL ,L 

'-U= / [u 1 (x) - u(x)]dx = I 

Jo JO 



(28) 



where 6u(x) and 8U represent the corresponding differences between the values of 
internal energy for the two set-ups. It is easy to realize that 8U and 6u(x] can only be zero 
when, 

u'(x(T')) = u(x(T}} (29) 

which can be true only if T(x) = T'(x). Alternatively, if c v of the bar is assumed as fairly 
constant then eq. (28) reads as 

f L 
8U = c v / 6T(x)6x, (30) 



where 8T (x} = T'(x] T(x}. Equation (30) clearly shows that 6U becomes zero only 
when T'(x] = T(x). 

Moreover, the preceding deductions remain true irrespective of V7" = T[ T[ = T\ 
T r = V7 1 or not. In the case of VT ^ VT" in general we have unequal values of the steady 
heat flux. Since in the case of a steady heat conduction q becomes position independent 
(that is throughout the length of the bar we have the same value of q) but since T remains 
position dependent makes u too a position dependent quantity. The above discussion 
clearly demonstrates that u is determined only by T and not by q. The latter quantity 
influences u via div q, that is by the net retention or rejection of energy. In the case of a 
finite but rapid change of 7s to T's and before a new steady state is reached we have 
T(x,t), VT(x,t) and q(jc, t). During this spontaneous attainment of the corresponding 
new steady state we still have u(x, t) = u(T(x, t}) and (du/dt)(x, t} is determined by the 
instantaneous value of divq(x, t) according to eq. (3) even if q(x, t) might be deter- 
mined by say MCV equation, eq. (1). When MCV equation operates there are two 



du_ d_d* FT 
p dt~ dx dt + *dx>' 

Notice that for the sake of simplicity r and A have been assumed constant of the operator 
<9/<9x. Further, on once more using eq. (3), eq. (31) gives 

d 2 u du . 8 2 T . 

rp W +p lTt- X M =Q - 

Notice that eq. (32) clearly shows the T dependence of u. Of course, on substituting (17) 
into (32) produces the well known hyperbolic equation of the telegraphist type, namely, 



where x = ^J pc v is the thermal diffusivity. This derivation of eq. (33) uses (17) in the last 
step and hence in this respect this derivation is different than the one generally followed 
in the literature in which eq. (32) gets skipped over. Recall that in deriving eq. (33) one 
uses (17) which is a companion to the Gibbs relation, eq. (16), and not to the extended 
Gibbs relation, eq. (10). That is eq. (33) is a description within the traditional Gibbs 
relation domain and not in the domain of an extended Gibbs relation. Hence we conclude 
that eq. (33) is not connected at all with the proposed extended Gibbs relation, eq. (21) of 
BIT. 

5. Systems in nonequilibrium and temperature 

This discussion will be incomplete if we do not clarify that what one is describing by T 
when a system is not in equilibrium. In BIT [1-7] in particular and in thermodynamic 
discussions [17,24,25] in general, there exist conflicting claims about the concept of 
temperature for systems in nonequilibrium states. The views expressed are based on both 
macroscopic and microscopic considerations. However, herein we are restricting 
ourselves to the macroscopic considerations only because thermodynamics operates at 
this level. As said above of eq. (21) is not believed in BIT as the same physical entity as 
T for an equilibrium case. For example from eq. (21) a so-called thermodynamic defini- 
tion of nonequilibrium temperature follows as d~~ l = (dr]/du) qv . At other places the 
symbol 9 is used and is referred to as the contact temperature [2427] and is described as 
the temperature of that heat reservoir which when brought in contact with the local 
nonequilibrium tiny system does not produce a heat flux. In view of these two divergent 
stands it becomes expedient to reexamine and understand the very basis of the 
experimentally measured temperatures. 

Recall that when one measures temperature it is required in principle that this act 
should not significantly alter the state of the system (or precisely that of the local pocket 
of the system in contact with the measuring device). For example, when a thermometer 
is inserted in a system in equilibrium the immediate reaction is that a finite heat flux is 



The quantity of energy utilized by the thermometer is insignificantly small so that the 
state of the system remains practically unaltered. Now carefully scrutinize this 
experiment and enquire when a steady reading on the thermometer scale will be 
registered. The answer is when the thermometric fluid stops the net acceptance (or 
rejection) of the energy from its immediate surroundings. That is when the internal 
energy of the thermometric fluid becomes steady. From eq. (3) the sufficient and 
necessary requirement for this is that the divergence of the heat flux within the 
thermometric fluid vanishes. Notice that the vanishing of the divergence of the heat flux 
does not demand the vanishing of the heat flux itself. However, it is the consequence of 
the equilibrium state of the system, which had no heat flux before the insertion of the 
thermometer, that simultaneously the heat flux within the thermometric fluid dies out. 

Now consider a system in nonequilibrium which at a given time and position has a 
finite value of heat flux. In this case too when a thermometer is inserted at a given 
position of the system the value of the heat flux at that position gets immediately altered 
by a significant amount. But after the passage of a very short period of time the 
thermometer shows a steady reading and the said position of the system practically 
regains its initial state. This happens due to the vanishing of the divergence of the heat 
flux within the thermometric fluid and not because of the dying out of the heat flux itself. 
If the latter would have been the case it would be at the expense of altering the state of the 
system. Thus we conclude that the primary requirement which is fulfilled by the 
temperature measurements is the vanishing of the divergence of the heat flux within the 
thermometric fluid. Hence, whether a system is in equilibrium or not the temperature 
measurements give us the identically same physical entity. Therefore, there is no such 
thing like nonequilibrium temperature (the adjective nonequilibrium applies only to a 
system and not to its state parameters) neither one is required to evoke a concept of 
contact temperature with the help of a hypothetical set of heat reservoirs. Since the 
definition of the contact temperature involves the condition of no heat flux then this 
cannot be achieved without altering the internal energy at the given nonequilibrium 
position of a system. Thus the contact temperature so obtained corresponds to a different 
state of the system and hence it is no wonder that this so-called contact temperature is 
obtained as independent of the internal energy of the tiny pocket under investigation 
[24,26]. That is, the contact temperature is not related at all to the given state of the 
system and hence is a hypothetical entity irrelavant for thermodynamic descriptions. It 
may be useful for a constitutive theory and hence for dynamics but this fact, in any way, 
should not dictate thermodynamics and its basic framework which is governed solely by 
the laws of thermodynamics. 

From the above discussion we establish that there is no distinction between the 
physical meaning of the temperature of an equilibrium and that of a nonequilibrium state 
of a system because the mechanism of registering temperature does not differ in the two 
cases. Indeed in a nonequilibrium state of a system there is a spatial temperature field 
while in the case of equilibrium the temperature is uniform throughout. The next obvious 
query is about the shortness of the time required for the dying out of the divergence of the 
heat flux within the thermometric fluid. This is determined by the corresponding 
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temperatures in a wide variety of nonequilibrium situations right from the industrial set- 
ups, laboratory experimentations to the biological systems. Hence to say that there is no 
well defined temperature in nonequilibrium tantamounts to shun our minds from the stark 
realities. 



6. Concluding remarks 

From the above discussion we see that q cannot assume the role of the a variable in eq. 
(14) and hence s cannot have q as an independent variable. Thus we once more conclude 
that the physical fluxes such as q etc. or the physical forces VT etc. cannot serve as 
independent thermodynamic variables to construct an extended thermodynamic frame- 
work for describing certain irreversible processes which are (believed) not covered by the 
traditional Gibbs relation, namely, 

^! = l^j_^^_ V^~ (34) 

dt~ Tdt Tdt ^ T dt ' ( ' 

A, 

The additional thermodynamic variables, as said above in 3, must find their origin in 
the additional nonequilibrium interactions not related with the three fundamental 
interactions namely thermal, mechanical (which includes both the short and long range 
interactions) and which cause the changes in chemical composition. The same 
conclusion, that is the heat flux etc. are not suitable to act as additional independent 
thermodynamic variables, has been earlier arrived at by us [20] on considering the 
dictates of the Gibbs-Duhem equation corresponding to eq. (6), the extended Gibbs 
relation of GPITT and also on the basis of the requirements of the Pfaffian exact 
differential form for the entropy function and by Garcfa-Colin et al [23] on the basis of 
kinetic theory. However, if one uses q etc. as independent variables then it amounts to 
tampering with the concept of entropy based on thermodynamic laws. That is then one is 
not talking about the Clausius entropy and hence not working within the domain of 
thermodynamics. 

Also we have seen above that an irreversible thermodynamic description contains in it 
the local thermometric temperature corresponding to a given nonequilibrium situation 
and none else. Of course, one is free to coincide it with the temperature on the Kelvin 
scale as one does in the case of equilibrium. The mechanism of registering temperatures 
involves the vanishing of the divergence of the heat flux within the tiny volume of the 
measuring device in contact with the local tiny pocket of the system and hence one 
measures the identically same entity as temperature in equilibrium as well as in non- 
equilibrium situations. 
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Abstract. We derive the metric for a Schwarzschild black hole with global monopole charge by 
relaxing asymptotic flatness of the Schwarzschild field. We then study the effect of global 
monopole charge on particle orbits and the Hawking radiation. It turns out that existence, 
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1. Introduction 

The Schwarzschild solution is the unique spherically symmetric solution of the vacuum 
Einstein equation. It represents a static black hole. It turns out [1] that a global monopole 
charge accompanies spontaneous breaking of global 0(3) symmetry into f/(l) in phase 
transitions in the Universe. Putting on global charge to the Schwarzschild black hole will 
amount to breaking the vacuum and asymptotic flatness of the Schwarzschild spacetime. 
It still represents a localized object with a horizon. Our main aim in this paper is to study 
the effect of global monopole charge on particle orbits and the Hawking radiation. 

The spacetime of pure global monopole charge, when the Schwarzschild mass vanishes, 
can be regarded in some sense to be 'minimally' curved for gravitational charge density, 
47rp c = RikU l u k , u'ui = 1, vanishes. That means the relativistic (active gravitational) mass 
is zero for the spacetime [2], All radial trajectories will always remain parallel. Thus 
introduction of plohal mononole charge does not significantlv alter the nature of the 



z. me metric 

Let us write the general spherically symmetric metric 

ds 2 = Bdt 2 - Adr 2 - r 2 (d0 2 + sin 2 Odip 2 ) , (2.1) 

where A and B are functions of r and t. 

Now the equations RQ\ = and R = R\ will imply B = A~ l = 1+20, where (f) = <f>(r) 
[4]. Note that no boundary condition has been used. With this 

R Q =R\=-V^ (2.2) 

and 

*2=*i = -^W, (2-3) 

where a dash denotes derivative with respect to r. Then RQ = gives the good old 
Laplace equation which has the well-known general solution 

= fc- . (2.4) 

Now R\Q will determine k = and we obtain the Schwarzschild solution. Retaining k 
will make R\^Q anc ^ w ^ gi ye rise to stresses, 

(2 ' 5) 



These are precisely the stresses at large r required for a global monopole as we shall see 
in the next section. Here k will be related to the global monopole charge which produces 
neither acceleration nor tidal acceleration for radially moving free particles. Note that the 
gravitational charge density, 4-7rp c = Riku l u k , continues to remain zero indicated by 
vanishing of RQQ. 

3. The gravitational field of a global monopole 

Several authors [1, 5, 6, and the references therein] have discussed the gravitational field 
of a global monopole formed by spontaneous symmetry breaking of a triplet of scalar 
fields with a global symmetry group 0(3). The Lagrangian density of the isoscalar triplet 
iff with a = 1,2,3 is 

^W fl ) 2 ~(W-7? 2 ) 2 . (3-1) 

Topologically non-trivial self-supporting solutions to this system can be found. The 
ansatz describing a monopole is 

-, (3.2) 



where 

*V = r 2 , (3.3) 

with .^-meaning the corresponding cartesian component of x. Outside the monopole core 
we must have 

/(r) -> 1 as |x| -* oo (3.4) 

for V(ip) to vanish asymptotically. The vanishing of energy density is however not fast 
enough because 

. (3-5) 



The stress tensor of the system outside the core can be approximated as 

tf = rl= (3.6) 

and the other components vanish. These are precisely the stresses generated (eq. (2.5)) by 
keeping k in (2.4). 

Thus the most general metric for 7$ = T\ ^ and the remaining stress components 
being zero is 



Ik- d/ 2 - + 2* - ~ -r 2 d0 2 -r 2 sin 2 0d^ 2 , (3.7) 

r J \ r ) 

with k = 47r?7 2 and M having the usual meaning of mass of the central gravitating object 
[1]. Here then central mass is that of the monopole which is usually negligible. However 
if we consider the spherically symmetric gravitational collapse of the matter around such 
a monopole, we find that a black hole is formed. The parameter M then corresponds to the 
mass of the black hole. 
Even when M = 0, this spacetime has non-zero curvature: 

RQ= R\ = ROI = (3.8) 

but 

*t = ^. (3.9) 

r 2 - 

That the total energy of such solutions is divergent makes them unrealistic except 
perhaps as fleeting entities in the course of a cosmic phase transition [1]. The core mass M 
is nevertheless finite. It is interesting to note that retention of k in (2.4) amounts to giving 
up asymptotic flatness and that exactly generates the stresses of a global monopole. 

A remarkable conclusion that follows from this is that a black hole can possess 'hair' in 
the form of topological charge. How will this charge be detected? Topologically 
nontrivial solutions are known to have the peculiar property that their symmetry 
generators are linear combinations of spacetime symmetries and internal symmetries. In 
the present case, it would mean that a motion on a sphere of constant radius r has to be 



Similar conclusions should follow for topological sectors of higher winding number, 
except that the spacetime metric will not be in the simple form corresponding to spherical 
symmetry. 

Another interesting manifestation of the topological charge can be obstruction to 
global definition of the wave function of certain particle species. For example, if the 
symmetry group is really the spin group SU(2} instead of 0(3} (more generally, spin (N} 
instead of O(N}} then an isospinor species x ma Y exist. Now a 2?r space rotation will 
induce a 2^0(3} rotation which will change the sign of the wave function which therefore 
cannot be globally defined. 

Note that the fate of the global charge S z of the 0(3} group will be the same as that of 
any other global charge such as the baryon number. 

4. Particle orbits 

In this section we wish to investigate the effect of global charge on the particle orbits; 
existence, boundedness and stability of circular orbits, light deflection and perihelion- 
shift. 

4.1 Circular orbits 

For the metric (3.7), we have 

At = E, r 2 p = l, (4.1) 

where A = 1 + 2k - 2M/r, t = dt/ds; E and / are specific energy and angular momentum 
of a test particle. Substituting (4.1) in the metric (we have set = n/2 as usual), 

i "> / A 1 i 

r 2 = 2 -A (- + l)=E 2 -V 2 . (4.2} 



For existence of circular orbit, we should have both r = r = 0. This will give the 
existence threshold 



where equality refers to the photon orbit. Of course / > 2\/3M/(l + 2k}. For bounded- 
ness, we should have E 2 > (1 + 2k}, which leads to 

" 



while stability of circular orbits will further require V" > 0, giving 
6M 



(4 ' 5) 



Following the standard calculation [7], we write 



-1/2 



dr 
r 



(4.6) 



where ro is the radius vector of the closest approach to the gravitating body. The 
deflection angle will be given by 

7T 



-ip(oo)\ - 



T+2F 



(4.7) 



Note that when M = 0, ip has the period 27r/vT+~2& and not 2-Tr, there is a deficit angle 
similar to the case of cosmic string. From (4.6) and (4.7) we get 



r (l+2fc) 32 H) 

4.3 Perihelion-shift 

Here the corresponding relation [7] is 



2 



^ 



A- 1 / 2 

x dr, 

r 2 

where r refer to aphelion and perihelion points. The perihelion shift is given by 



(4.9) 



T+2& 



and we thus obtain 



L(l+2ky/2 



(4.10) 



(4.11) 



where 



2 



Thus effect of global charge is to scale the Schwarzschild values by (1 STTTj 2 )" 1 for 
existence, boundedness and stability thresholds of circular orbits and by (1 &7rr) 2 }~ 3 ' 2 
for light deflection and perihelion-shift. It means a slight enhancement in the 
Schwarzschild values. 

5. Semiclassical effects 

Next we investigate the question of Hawking radiation in this background metric. Since 



Yi m is a spherical harmonic and R^i satisfies the equation 



Consider the equation satisfied by a massless scalar quantum field i/j 

n^ = o. (5.1) 

The advanced and retarded coordinates for this spacetime are 

= f + r', (5.2) 

u = t-r\ (5.3) 

where r* is the tortoise coordinate defined in this case as 

, rdr . N 

dr = . (5-4) 

1 + 2k (2M/r) 

The wave equation can be separated in (?, r, 9, </?) coordinates. The mode functions 
relevant to this spacetime are of the form 

D /,.N 

(5.5) 

(-Jj-*2 ' V l*^* ' *~J' ' ~" p *' JL i ' "' v <-"" \) M V / 

Because of the 'potential' term in square brackets, the standard incoming waves will 
partially scatter back the gravitational field to become a superposition of incoming and 
outgoing waves. However, the effective potential vanishes as r* * 00 (i.e. for 
r > 2M, oo). Thus, in those regions the mode functions will be of the form 

and fie-^yA (5.7) 

The surface gravity of the black hole formed as a result of the gravitational collapse is 

^(1+2/0^(1-87^ ( ^ g) 

This redefinition of the surface gravity reduces the calculation to that of the 
Schwarzschild case. With this redefinition, the Penrose diagram of the present system 
is identical to that of the Schwarzschild case. 

As in the Schwarzschild case, assuming there are no bound states at late times, the 
ingoing and outgoing modes form a complete basis in terms of which the field at late 
times can be expanded. Then an outgoing mode of the form r~ 1 e~' WM y/ m at late times can 
be traced back to early times in the same fashion as in the Schwarzschild case. At early 
times, the ray, moving along constant phase v lines will have the form e~ iu}U (v). The affine 
parameter on the past horizon can be chosen as 

A=-ce . (5.9) 

Then the function 



(5.11) 

V ' 



Thus the thermal spectrum is recovered but the temperature is now (1 - 87n7 2 ) 2 /87rM. 

Consider next the quanta of the ijj field responsible for the monopole. The equations 
satisfied by 

i> a =^ a -'tt g (5.13) 

are, after linearization, 




= 0. (5.14) 

Asymptotically the equation becomes 

= 0. (5.15) 




Consider the ^=3 mode and the asymptotic z direction, the other two (x and v) remaining 
finite. Then the equations become 



+ IXrf = 0, (5.16) 

=0, (5.17) 

= 0. (5.18) 

This corresponds to the radial mode being massive but the two transverse modes being 
massless. Thus the preceding analysis will apply to the two transverse modes as well. We 
see that an imprint of the topological charge of the black hole is left on the Hawking 
radiation as well, namely under angular displacement, different components of the t/> 
triplet will be detected. 

Thus the fate of our black hole is similar to that of the magnetically charged black 
holes [8]. While the question of its ultimate fate cannot be settled within the semiclassical 
approximation, the answer will be the same as for a Schwarzschild black hole. If it can in 
fact evaporate completely, and does so, then it leaves behind a monopole identical to the 
original one before collapse. 

In concluding this section, we note that the persistence of the thermal radiation is not 
surprising. It is known from the work in axiomatic field theory [9] that when horizons 
exist, a unique nontrivial automorphism exists for the ground state. The automorphism 
works out to be periodicity in time, leading to the Kubo-Martin-Schwinger boundary 
condition on the Green function. There has also been extensive work on thermodynamics 
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6. Discussion 

The main purpose of this investigation was to study the effect of global monopole charge 
on particle orbits and the Hawking radiation. It turns out that existence, boundedness and 
stability threshold for circular orbits scale up by (1 STrr? 2 )" 1 , perihelion shift and light 
bending by (1 STrr/ 2 )" 3 , while the Hawking temperature scales down by (1 STrry 2 ) 
the Schwarzschild values. This is how the global monopole charge r\ affects the particle 
orbits and the Hawking radiation thermal spectrum. It may be noted that the event horizon 
is given by r = 2M(1 STr?? 2 )" 1 , while the red-shifted proper acceleration remains 
unchanged as M/r 2 . It is the scaling up of the horizon that leads to scaling down of the 
surface gravity (temperature) of the hole. 
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Note added in the proof: 

By resolving the Riemann curvature into electric and magnetic parts, it has now been 
shown [12] the Schwarzschild black hole with global monopole (3.7) is dual to the 
Schwarzschild black hole without global monopole. The solution (3.7) with M = is 
dual to flat spacetime and hence it can justifiably be considered as "minimally" curved. 
By duality we mean the interchange between active and passive electric parts of the field. 
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In the small -x limit, the coupled DGLAP [1,2] sea-quark and gluon evolution equations 
for leading order can be written as sea-quark and gluon evolution equations in decoupled 
forms. The gluon evolution equation is of the form of wave equation in terms of the 
variables (= In(jc /*)) and (= ln[ln(g 2 /A 2 )/ln(g/A 2 )]) (X Q and go being the starting 
scales of perturbative evolution) [3]. In the asymptotic limit of a(= -\/|C), the solution of 
this wave equation gives rise to scaling in terms of a and./?(= \//C) which is termed as 
'double asymptotic scaling' (DAS) [3]. From the gluon momentum distribution function 
we can find the sea-quark momentum distribution function using DGLAP sea-quark 
evolution equation for small-jc. Thus, we can obtain the structure function. The prediction 
for proton structure function from DAS, when compared with the data from HERA [7], 
has good agreement [3-5] for large values of cr and p. In the case of gluon, prediction 
from DAS agrees well with the fit obtained by Gliick et al [6] above a certain value of 
p(= 8/y}, which can be found analytically [8]. Hence, we found a simple solution [9] of 
gluon wave equation at fixed p which agrees well with the fit [6] below p = 8/y. Here, we 
extend the solution [9] in approximate form to proton structure function for large value of 
C (that is, Q 2 ) which can be used to analyse the data below the critical value of p(= 8/y}, 
because, the complete expression for structure function from DAS has a discontinuity at 
p = 8/y and is negative below this value. To compare our prediction with that from DAS, 
we separated the HERA data [11] into two sets one set having p > 8/y and another one 
having p < 8/y. The DAS predictions, without inclusion of complete splitting kernel of 
the process g > qq, for proton structure function is found to be in good agreement with 

rloto Knt nut- rvr<a/-lirti<-\o fall -For Halrtti/ tV>o rlata fnr n ^* A Ai T5nt frvr n ^" X I <\\ r\v\r 



where G(, C) = xg(x : Q 2 } is the gluon momentum distribution function, y = y/12/A), 
6= (11+ 2n//27)/j0b, /3b = 1 1 - 2n//3, and n f is number of flavours. 
When p is fixed, the solution of (1) for large is [9] 



G(Cp 2 , C) - ^ exp[M + v^+WX/2]. (2) 

It is well known that, in the small-* limit, the contribution from gluon is very large 
compared to quark contribution. Hence the LO-DGLAP evolution equation for sea-quark 
can be approximated as 



- 47T 
or 



(3) 

where q(x, Q 2 ) is the quark momentum distribution function, P qg is the leading order 
splitting kernel for g * qq process and '&>' indicates convolution with respect to the 
variable xf. a s is the strong coupling constant. 
Since, in the LO 

^Hi[z 2 +(i-z) 2 ] (4) 

hence, (3) can be written as 



tt, C) = dt>\2*-Q + 1 - 2e]G(, C). (5) 

^C 6 Jo 

In the small-* h'mit, the proton structure function can be approximated as 

*S&C)=tk$&0 (6) 

Therefore, using (2) in (5), the expression for F% can be written as 

p _5A_ 2 J 

F2 ~ n/7 P 



-fa - 2) (p 2 



2a 
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-ru)\ir -i- a - i; p~(p-fa)j 



where a = [(-5 + ^/<5 2 + 4j> 2 p 2 )/2]. 
When ~ 1 



108 



2a 2 la a(a + 2) 2(a + 2) a(\ + a) l+a 
If is large, 

-, P / j. 9 j-\ J-^ "7 nf I A 1 -I 



_ 5A 2 g< I" a 2 + a + 2 
~ 108 n/J T L2a(a+l)(a + 

For fixed p, a is fixed and hence 

*?="y 2 ^ (7) 

where 

5A 



' 108 ^(a + !)(* + 2)J' 

The normalising factor N has a dependency on p. But, it has no singularity for all p > 0, 
which is the advantage over the DAS result [3]. 
From (7) we can have 

f?~tfj> (8) 

the approximate double scaling, when p is large, which is a consequence of double log 
approximation [10] (it is due to gluon density). For fixed p, the slope of the curve, In 
F% vs. cr, is half to that predicted by DAS. It shows that for approximate double scaling 
of F it is not necessary to use DAS form of splitting kernel for the evolution of sea- 
quark. In both the cases subasymptotic contribution from quark (P qq ] is neglected. 

To compare with HERA data [11], we use (7). In deriving the expression for proton 
structure function, our approximation is such that it is applicable for large , that is, 
large Q 2 .. 

It has been observed in gluon that [8], there is a lower bound on p, below which the 
prediction from DAS differs from fit obtained by Gliick et al [6]. To find a similar 
behaviour in proton structure function, we use the complete form of F%, that is 



C is a constant, y = y/p, G(cr, p) = cr 1/2 exp I 2^a 8- j , 



I. /,00= 1--5T 



-i 



T>V 2(l+y)(2+y)(3+y) 



and the limiting form of (9) (when 3; > 0), 
Ff(cr, p) - C -a 1 / 2 exp /2yer - 



(10) 



<T, p, ^ and C are evaluated using XQ = 0.1, <2o = 1 GeV 2 A = 0.23 GeV. Ff is evaluated 
using nf = 4 and = 5. 

It can be observed (figure 1) that eq. (9) has a singularity at p = <5/y(~ 1.13) due to 
f q (y) and uncertainty of this function is much greater than/ g (y) near the 'singularity' and 
hence its effect cannot be neglected. It is easy to find that F% is negative according to eq. 
(9) for p < 8/y. The prediction for F\ from DAS in the region p < 8/y contradicts QCD 
as in the case of gluon [8]. So, p = 8/y is being termed as 'critical point'. 

Next, we separated the HERA data [11] into two sets (figures 2 and 3). One set is for 
p < 1.13 (indicated by traingles) and another one for p > 1.13 (indicated by circles). The 
shaded regions (with dots) indicate the average error in experimental data. In figure 2 the 
dotted-dashed line indicates the prediction from DAS for p = 2 (using (10)) and the solid 
line is the prediction from (7) for the same value of p. For the predictions, we put 
A -1.4763 and C = 0.1245. The prediction from DAS agrees with data. But, the 
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Figure 2. In F P 2 plotted against cr using eq. (10) (the dotted dashed line) and eq. (7) 
(the solid line) for p = 2. The black circles indicate experimental data points having 
p > 1.13 and the triangles indicate data points having p < 1.13. The region shaded 
with dots indicates average uncertainty in experimental HERA data. 
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Figure 3. In F plotted against cr using eq. (7) (the solid line) and eq. (10) (the 
dotted dashed line) for p=l. The black triangles indicate experimental data points 
having p < 1.13 and circles indicate data points having p > 1.13. The region shaded 



found that for p < \.\3,f q (y] ^ 1 and so, use of (10) in this region is questionable. 




Figure 4. Ff vs. jc for <2 2 = 20,25,35, 120,200,250,350 and 500 GeV 2 . The 
predictions from eq. (9) (the solid lines), eq. (10) (the dotted lines) and eq. (7) (the 
dashed lines). Data points with error bars are from HERA [11]. 



such that variation of p is very small, which may be considered as fixed. Under this 
consideration (7) can be used for conventional plot of x vs. F at fixed, but large Q 2 . In 
figure 4, we plot x vs. F P 2 from (7) (dashed line), (9) (solid line) and (10) (dotted line) for 
Q 2 = 20, 25, 35, 120, 200, 250, 350 and 500 GeV 2 . Equation (9) agrees well with data for 
Q 2 = 20, 25 and 35 GeV 2 (figure 4 (i,ii,iii)). But for large Q 2 (Q 2 > 120 GeV 2 ) it does not 
agree with data for comparatively large values of* (figure 4 (iv,v,vi,vii,viii)). Agreement 
with data is found for 56 data points out of 169. Equation (10) agrees well with data for 
all the Q 2 and nearly all x (agreement is found for 1 15 points out of 169). Equation (7) is 
found to be valid for Q 2 = 120,200,250,350 and 500 GeV 2 and x > 0.02 (figure 4 
(iv,v,vi,vii,viii)) (agreement is found for 45 points out of 169). Keeping p ~ 1 as we go on 
increasing Q 2 ,x can be decreased and (7) can be used to smaller values of x. 

Thus, we have shown that DAS is not the only one solution to explore the physics at 
HERA. There is a critical value of p below which application of DAS is questionable. The 
difference in slope of the curve InFjVS.cr for data having p < 8/y with the data having 
p > 8/y indicates a change in physics of low-x, that is, a change from hard pomeron to 
soft pomeron [9]. More data below and above the critical value of p will properly indicate 
the change. 
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Abstract. The weak radiative decay A/, > A^ is studied in the heavy quark effective theory 
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1. Introduction 

The weak radiative decays (flavor changing radiative decays) of hyperons have proven to 
be a challenge to both theorists and experimentalists. Experimental difficulties result 
from small branching ratios of these decays and their copious photon background. These 
transitions require the joint occurence of weak and electromagnetic interactions. Their 
theoretical treatment requires the inclusion of separate short distance (SD) and long 
distance (LD) contributions [1,2]. The relative size of the two types of processes is an 
issue to be determined for every specific process. At the quark level there are three types 
of processes which contribute to the weak radiative decays of baryons classified as single, 
two and three quark transitions. The two quark transition corresponds to W exchange with 
the photon radiated by the participating quarks. The three quark transition, where the 
quark not participating in W exchange, radiates the photon, which is strongly suppressed. 
The single quark transition involves a (SD) contribution due to electromagnetic penguin 
diagrams [2, 3]. The weak radiative decay A fo * A.J is described by the single quark 
transition and is dominated by the short distance electromagnetic penguin transition 
b > sj. The recent measurement by CLEO of B > AT* 7 [4] and B > Xj [5] processes, 
confirm that this mode is indeed dominated by the electromagnetic penguins. 

After knowing the fact that the decay process A& > A7 is dominated by the short 
distance electromagnetic penguin, we use heavy quark effective theory to study this decay 
Drocess. Important nroeress in the theoretical description of hadrons containing one 



particles involved, and is normalized at the point of zero recoil. Similarly in case of weak 
decays of heavy baryons one can also write the form factors in terms of another Isgur- 
Wise function [8, 10]. 

The idea of heavy quark symmetry is considered to be effective only for heavy quarks 
whose masses mQ are significantly larger than the QCD scale AQCD ( 200 MeV). Thus, 
the hadrons including c/b quarks are reasonable laboratory to test the heavy quark 
effective theory. The s quark is not heavy but not light as the u and d quarks. In many 
cases the light flavor 517(3) symmetry is successful for s quark system after 
nondegeneracy of .s-quark mass is taken into account as symmetry breaking. However 
some features of s quark physics are understood if the 5 quark is treated as heavy. A 
notable example is the Okubo-Zweig-Iizuka suppression of the decay (f> TTTTTT, which is 
a direct consequence of the fact that the s quark mass or twice the 5 quark mass is heavier 
than the QCD scale. The s-quark seems to be marginal in both heavy quark symmetry and 
light flavor one. Therefore it is interesting to examine the applicability of the heavy flavor 
symmetry with corrections, for the .s-quark system even if s-quark mass is not large 
enough compared to AQCD- Although this assumption might be too optimistic the previous 
challenging studies [1 1-14] suggest that heavy quark effective theory is still effective for 
s-quark system if the constituent quark mass is employed for m s and l/m s corrections are 
taken into account properly. The 1/mq corrections in the HQET is formulated by putting 
the heavy quark Q almost on on-shell and expanding in powers of the residual momentum 
k which is generally written as AQCD/^Q- Then with AQCD = 200 MeV and m s = 550 MeV, 
a measure for l/m s correction becomes AQCD/WJ 0.36, which is expected to keep the 
perturbation meaningful. Treating the s quark as heavy Ali et al [13] have studied the 
radiative rare decays B > K*^. Ahmady and Liu [14] estimated the B K^f decay rates 
and obtained quite good results. However in a recent letter Chakraveity et al [15] have 
shown that the application of heavy quark symmetry is not reliable for .s-quark system. 
On the other hand applying HQET to B > K^ transitions, Roberts and Ledroit [16] have 
suggested that the heavy s-limit may give an acceptable description of unpolarized data, 
but does not give reliable results when applied to polarization observables. 

In this work we intend to test the applicability of heavy quark effective theory to the 
weak radiative decay Aj, A.J including corrections of the order 1/mQ. In the HQET there 
are two types of 1/mQ corrections to the hadronic form factors: one from 1/m.q corrections 
to the current operators and the other from the presence of higher dimensional operators in 
the effective Lagrangian [9,17], The latter amounts to the hadronic wave function 
modification. As shown in [17] all the form factors at 1/mQ order are still expressible in 
terms of a single function and one dimensional constant. Thus including 1/mQ corrections 
arising from both the sources, we calculate the branching ratio for the decay process. The 
Isgur-Wise function is evaluated in the large N c limit [18]. This decay is also studied in 
[19] in the effective Lagrangian approach, where they have incorporated two ideas - the 
heavy quark symmetry scheme with both b and s treated as heavy and the MIT bag model. 
However, regarding the 1/mQ corrections, they have included the corrections of order 
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our predicted result should be much more reliable than that of theirs. 

The paper is organized as follows. In 2 we present the general framework for the 
study of the radiative rare decays. The Isgur-Wise function is evaluated in the large N c 
limit in 3. Section 4 contains results and discussion. 

2. General framework 

The short distance electromagnetic penguin process is governed by the electromagnetic 
penguins operators. At the weak scale ^ = m b , the transition b > 57 is mediated by the 
effective Hamiltonian [20] 

, fe v,;c 7 K)o 7 (m,), (i) 

where the operator O 7 (mi) is given by 

v (l - j 5 )b] (2) 



and the Wilson coefficient C 7 (m&) is 

C 7 (m 4 ) = rT 16723 ^^) + tl(>7 10/23 - 1) +i>(^ 8/23 - 1)1 (3) 

with 77 = a s (m b }/a s (M w ) and 

x(to?+5x-1 3^-2jc 7 \ m? //n 

-^ x= - (4) 



The value of C 7 (mfc) is taken as C 7 (m^)= 0.31 [21]: The transition amplitude for the 
decay process A^, > A7 is given as 



M(K b -* A 7 ) = (A 7 |Hef f |A,) - ^y f ;c 7 (m i )(A 7 |0 7 |A fe ). (5) 

V2 

We shall use HQET to compute this matrix element. In the HQET the heavy A Q baryons 
are described as the bound state of the heavy quark field and the light diquark field. In 
terms of the interpolating fields it is given as [22] 

A Q (v, 5) = MA Q (v, s)(j> v h v , (6) 

where c/> v represents the light diquark field and h v is the heavy quark field. MA Q (V, s) is the 
usual Dirac spinor, which obeys the constraint relation 

j/tu(v,s} = u(v,s}. (7) 

The heavy quark field in the effective theory hq(v,x) is related to the quark field Q(x) of 
the QCD Lagrangian by 



Q'rQ = h'Th + - hT(ih + AHW^, (9) 

2/7Z 



where h and h' are the heavy quark field operators in the effective theory for the quarks Q 
and Q'. In the effective theory the matrix elements of the operators on the right hand side 
of (9) between the Aq-baryon states are given as [23] 

(10) 

, (ii) 

,.0. (12) 
The most general decomposition of a involves two scalar functions defined by 

( a (v,v'}=( + ( W }(v + v') a + ^( W )(v-V f ) a , (13) 

where w = v-if, C(w) are given as [23] 



C-H=CH- (is) 

The parameter A is the scale of light degrees of freedom given as A = mA Q - mq 
700MeV. Now at 1/mq order arising from the current operator modification, the 
matrix element of the operator O? between the A;, baryon in the initial state and a photon 
with A baryon in the final state may be written as 



l ( 7 < "A i ' v ^ \ i *A 1 "' > ** 1 l u " l- 1 ^ 1,^5 ^// ) v ^/ 

where 

^_-w[ / ! l 



(17) 
and e^ is the polarization vector of the emitted photon. Using eqs (10-17) one can obtain 



/a(1 +75)7Q 
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2m s \w + i w + 1 



_ 75) 



. 



As for the corrections arising from the hadronic states, we must also include I/WQ 
corrections to the effective Lagrangian, i.e., 



where 

Now inserting the higher order term of the effective Lagrangian ({) into the matrix 
elements of the lowest order current J = h'Th, we obtain [23] 

I f ? ^\~\~^ ( \ /*0 1 \ 

Insertion of the kinetic operator O\ preserves the Dirac structure of the current and trans- 
forms trivially under the spin symmetry, contributing to the form factors in the same 
proportion as the leading order term (10) and hence effectively correcting the Isgur-Wise 
function [9, 17, 23]. The contributions from the chromomagnetic operator Oi vanishes by 
Lorentz invariance. For the subleading form factor, vector current conservation implies 

(22) 



Since x( w ] on ly corrects the Isgur-Wise function, we assume a simple form being 
consistent with the constraint (22) as done in [12] for the mesonic case 

. (23) 



Thus with eqs (5, 18-23) we obtain the transition amplitude for the weak radiative decay 
A A7, including corrections of order (!/WIQ) as 



= iu A (v',s'}(A + S 75 )< W ^'' MAb (M, (24) 
where A and B are given as 



(25) 

and' 

GP e . . . ........ T_ /A AM 



\kf(\A\ 2 + |5| 2 ), (27) 

where k is the c.o.m. momentum of the emitted photon. 

Since it is not possible to find out the Isgur-Wise function in HQET, we will evaluate 
the same using some model. Recently a breakthrough has been achieved to evaluate the 
IW function for baryons i.e., the bound state soliton picture in the large -N c limit [18]. 

3. Evaluation of the Isgur-Wise function 

Here we have presented the evaluation of the Isgur-Wise function in the same manner as 
suggested in [18]. In the large N c limit the light baryons n, p, A etc; can be viewed as 
solitons in the chiral Lagrangian for pion self interaction [24]. The baryons containing a 
single heavy strange (or bottom) quark are bound states of these solitons with K and K* 
(or B and 5*) mesons [25-29]. In this paper we use the bound state soliton picture 
to estimate the value of the baryonic Isgur-Wise function. In the ground state of 
AQ baryons, the light quarks are in the spin state [22]. Hence in the bound state soliton 
picture, Aq-type bound state arise when the spin of the light degrees of freedom of the 
heavy meson and the spin of the nucleon are combined into a spin zero configuration 
where as the isospin of the heavy meson and that of the nucleon are combined into an 
isospin zero state. Other baryons (e.g., the A) only contribute to the bound states with 
higher isospin. 

Let the light degrees of freedom of the heavy baryon be denoted by |/, /a; 57, m/), where 
7 and si denote their isospin and spin quantum numbers while /3 and m/ are their three 
components respectively. Hence the light degrees of freedom of AQ baryon is denoted by 
|0, 0;0, 0). The chiral soliton is denoted by \R, b\ 7?, ), where R = 1/2 for the nucleon. 
On the other hand the light degrees of freedom of the heavy meson is given as 



In the large N c limit, the binding potential between the chiral soliton and heavy meson 
is independent of both the isospin and spin of the particles. Hence for the light degrees of 
freedom of AQ baryon, we have the decomposition as 

|0, 0; 0, 0(t;)) = J d 3 q<&(q)(l/2, *; 1/2, c|0, 0)(l/2, n; l/2,p|0, 0) 

, (28) 



where (7i,mi;J2,m2|J r ,M)'s are the Clebsch-Gordan coefficients. $(q) is the ground 
state wave function and M B and M# are the masses of the chiral soliton and heavy meson 
respectively. 

The spin- 1/2 A Q baryon is composed of a spin- 1/2 heavy quark and spin-0 light degrees 
of freedom. Hence the matrix element of the current h^(v'}Th^(v} between A^ and A 
baryons is given as 

C^^^^OjO.O (t/)|0,0;0,0 (v}}u A (v f ,s')Tu^(v, S }. 



(t;-t/) = (0,0;0,0(t/)|0 J 0;0,0(i;)) 

= /Vq\l 3 q$*(qO*(q)(l/2,tf; 1/2X10,0)* 

/ 

x (1/2, n'; l/2,/'|0, 0)*(l/2, ft; 1/2, c|0, 0)(l/2, n; l/2,/>|0, 0) 



(30) 

Using the normalization conditions for the chiral soliton and heavy meson states, all the 
Clebsch-Gordan coefficients in (30) turn out to be unity and the Isgur-Wise function is 
given as 

C(w -v')= f d 3 q$*(q)*(q + M B (v - v')). (31) 

It is noted from eq. (31) that the IW function depends on the spatial wave function $(q) 
of the AQ baryon. In the large N c limit, the binding potential between the heavy meson 
and the chiral soliton is simple harmonic [29], and hence the wave function is taken as 

(32) 



K is the spring constant and its value is taken to be (440 Me V) [30]. In the rest frame of 
the initial state, v = (1,0) and v' directed along z-axis we obtain the Isgur-Wise function 
(31) using (32) for non-relativistic recoils i.e., v'| 2 w 2(v v' - 1), as 

.. / (vv'-l) /Ml\ ,, 

c(W) = P M = -v-f (33) 

V 2 V K J 

The product (v i/) is determined by considering the kinematics of the system. Since we 
are dealing with the two body decays A.b(v) > A(v'}^(k}, momemtum conservation 
implies 

+ F, (34) 



which gives 

vv ! = M },". + ^. (35) 



With eqs (33) and (35), the value of the Isgur-Wise function is extracted as 

((v v') = 0.0794. (36) 

4. Results and discussion 

To estimate the numerical result we have used the following values. The values of the 



B(A b - A 7 ) = 1.48 x 1(T 5 . (37) 

In this paper treating s-quark as heavy an attempt is made to study the rare radiative 
decay A^, A 7 in the heavy quark effective theory. The decay is described by the single 
quark transition and is dominated by the short distance electromagnetic penguins. In the 
effective theory this decay matrix element can be written by an expression involving the 
Isgur-Wise function. Since the quark masses are not infinitely heavy, we have considered 
l/mQ(Q = b,s] corrections arising from the modification of current as well as the 
hadronic states. As suggested in [17] all the hadronic form factors at. order (1/mq) are 
expressed in terms of a single function, the Isgur-Wise function. Since HQET does not 
predict the shape of the IW function, it makes sense to extract the same using a model. 
Here we have evaluated the same using the bound state soliton picture, in the large N c 
limit [18] and obtained the branching ratio for the decay Ab > Aj. The present status of 
the subject, in particular the still unavailable data, do not allow for a more complete 
theoretical approach, however we are hopeful that future experiments will lead to more 
theoretical understanding. The predicted branching ratio for this decay is (1 0.5) x 10~ 5 
[32]. In an earlier investigation Cheng et al [19] obtained this branching ratio using heavy 
quark symmetry, MIT bag model and l/m s correction only to the current operator, to be 
(0.8-1.5) xlO" 5 . Additional experimental data would greatly help for a better 
understanding of the weak radiative decays of heavy hyperons, which can serve as a 
signal for new physics beyond the standard model. Our result indicates that heavy quark 
symmetry does work in the first approximation to calculate the exclusive decays 
involving b * s transitions. 
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Abstract. We address the problem of classical frictional motion under a potential V possessing a 
barrier, apart from other possible confining and nonstationary terms. It is pointed out that the 
Green's solution of the exact equation of motion can be reduced (under suitable conditions) either to 
an improved Rayleigh form or a non-Rayleigh form, the latter being outside the scope of the 
standard large-friction treatment of the Fokker-Planck equation. The resulting dissipationiess 
dynamics involves an appropriately scaled potential which may have promising applications to 
quantum stochastic phenomena. Genuine dissipative corrections in regions far away from the 
barrier can be accounted for by the higher-order terms in our asymptotic expansions. 
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1. Introduction 

The properties of single-particle, nonconservative, frictional motion in classical [1], 
quantum [2], stochastic [3] and statistical [4] mechanics are very important conceptually 
and quite interesting application-wise. Damping effects on classical trajectories are 
theoretically introduced via a linear- velocity term [la] in an extension of Newton's law or 
a quadratic dissipation function [Ib] in a generalization of Lagrange's equations or an 
explicitly time-dependent factor [Ic] in the Bateman-Caldirola-kanai (BCK) Hamilto- 
nian. Numerous applications of frictional trajectories include the calculation of 
deterministic collision [Id] of nuclear heavy ions and probabilistic description of 
phenomena-like Brownian movement [3a] of the harmonic oscillator, signal-to-noise 
ratio for stochastic resonance [3b] in a double well etc. It may be noted that the problems 
tackled in refs, [Id] and [3b] involve a potential energy which has a pronounced barrier 
(see figure 1) apart from other possible nonlinear/nonstationary terms. 

In the above context, previous workers have employed two different types of theore- 
tical modelling, viz., the standard frictional equation and its overdamped Rayleigh 
version. We review their salient features in 2 and point out that generally these equa- 
tions have been solved [Id, 3b] by numerical computation/analog simulation and their 
associated Lagrangians are inconvenient to use in practice. In 3 we carefully set up the 
formal Green's solution of the exact equation of motion and show analytically that the 




Figure 1. Schematic plot of the potential energy as a function of position (at a 
specific time). The barrier is explicitly shown. The dashed portions represent other 
confining and nonstationary parts of the potential. 



non-Newtonian dissipational trajectory near the barrier approximately coincides with a 
Newtonian dissipationless trajectory provided the applied potential is scaled by a suitable 
factor. Our concluding remarks appear in 4 where it is emphasized that a Lagrangian 
employing such a scaled potential can be useful both for path-integral based quantization 
and for the treatment of quantum stochastic resonance. The algebraic details of our 
formulation are relegated to the Appendix for convenience. 

2. Review of existing models 



Consider a test particle moving in an environment and let the symbols 
m, t , jc, v = djc/df , V, F = -dV/dx, y 



(la) 



respectively denote the mass, time, position, velocity, external potential, applied force 
and the coefficient of friction. This V may contain possible nonlinear terms in position 
and/or nonstationary terms in time. The symbols 



V f = 



= -dVj(x}/dx 



(Ib) 



will refer to the particular situation when the potential and force do not contain the time 
explicitly. 



2.1 Standard factional equation 

The basic equation of motion for the unknown trajectory x(t] valid for arbitrary damping, 
reads 



dv/dt + yV = F/m 
subject to the initial conditions 



(2a) 



Item Standard frictional model 



Motivation Dissipative force proportional to velocity is motivated experimentally 

Parameters The barrier frequency fi and the coefficient of friction 7 are arbitrary 

Estimated velocity For general potentials with a barrier v is obtained by numerical 

computation/analog simulation 

dv F 

Effective acceleration = yv + 

dt m 

Branch of motion Both types of branches are present depending on whether the particle 

is sliding down (Rayleigh-like branch) or climbing up (non- 
Rayleigh-like branch) the barrier 

Lagrangian L ^'(mv 2 - V) due to Bateman-Caldirola-Kanai which is ex- 

plicitly time-dependent 

Mechanical energy E, = mv 2 /2 + V, dE,/dt = -myv 2 < 



As far as x as a function of t is concerned, eq. (2a) is a second-order differential equation 
of the deterministic (stochastic) type according as the temporal of the force F is definite 
(random). Equation (2a) is equivalent to the formal expression 

dt' e^f/m, (2c) 

o 

where F' = F(x(t"), t'} is the force at the integration time t' when the position becomes 
x(f}. Other known features of the standard frictional equation are summarized in table 1. 
It is seen that irreversible dissipation of the mechanical energy E/ = mv 2 /2 + V/ is an 
inherent property of the model so that a dissipationless approximation to the motion is 
never attempted. Furthermore, since the underlying BCK [lc,2a] Lagrangian 
L = (mv 2 /2 - V] e yt is explicitly time-dependent, its path integration in general, would 
require complicated numerical matrix multiplication [5a]. 

2.2 Overdamped Rayleigh equation 

If the coefficient of friction is large enough, the right-hand-side of eq. (2c) can be 
evaluated by repeated partial integration to yield 

v w F/my - F/my 2 H ; yt > 1 ( 3a ) 

where F = dF/dt. On the so-called Rayleigh branch of motion (labelled by the 
superscript R), the leading term of eq. (3a) is picked up, giving a velocity parallel to the 
applied force, 

(*> - F/my. ( 3b ) 

As regards the trajectory x^ as a function of t, this is a first-order differential equation of 
the deterministic (stochastic) variety depending on whether F has a definite (random) 



Item Conventional Rayleigh model (label R~) 

Motivation The frictional force balances the applied force 

Parameter 7 large, yt 1 

p 

Estimated velocity v^ = 

my 

Effective acceleration - ~ W R ' H 

dt m 

Branch of motion Emphasis is on those branches where v and F are almost parallel as 

happens when the particle is sliding down the barrier in figure 1 

Lagrangian L^ does not exist because the acceleration has become trivial 

Potential energy vf^fc) such that dV^ /dt = -myv^ 2 < 



equation. It is observed that there is no recipe available to handle the so-called non- 
Rayleigh branch on which the velocity would be antiparallel to the external force. 
Furthermore, eq. (3b) being of the first order can not follow from a usual Lagrangian in 
which x^ is the only degree of freedom. 



3. The dissipationless Newtonian model 

3.1 Preliminaries 

In view of what has been said above it is worth examining (for finite 7) questions such as 
the possibility of a dissipationless picture, existence of a non-Rayleigh branch, construc- 
tion of a convenient Lagrangian, etc. This task will be accomplished analytically for those 
problems [Id, 3b] in which the input potential includes a parabolic barrier. Then the full 
external force on the particle can be decomposed as 

F = F B + F p = mtfx + F p (x, t) (4) 

where FB is the linear force due to the barrier, f2 the corresponding angular frequency and 
all the remaining contributions are lumped together in the perturbation force F p which is 
supposed to be confining at large distances. In the sequel we shall also need the useful 
symbols 

(5) 



3.2 Formulation 

The standard frictional equation (2a) is rewritten as 



t ' u<u. HJLIO iccuio LW a pan ui CAa^u, luiuicu iciaiiuna llivuivuig IlUcgiillS ^Scc /\ppendlX, 

eqs (A8), (A9a, b) between i>, x and F p , which can be suitably approximated depending 
upon the mechanical branch of interest. 



3.3 The 'improved' Rayleigh branch 

Suppose the particle is sliding down the potential hill of figure 1 so that its velocity is 
essentially parallel to the applied force. The particle then takes relatively less time to 
cover a given distance down the hill. After the system has relaxed with respect to the time 
scale il~' one employs the asymptotic expansion (All) and finds for the velocity 



---- ; n+fl. (7a) 

Here the dots . . . stand for nonleading terms of order F p /mfi?j_ and the superscript (R)* 
labels the Rayleigh branch as improved by our analysis. Of course, in the limit y/fi > 1 
the above result coincides with the usual Rayleigh estimate F/my. The improved 
expression for the acceleration is given by 



m 



which indeed has the appearance of a dissipationless Newtonian equation under a force 
that has been scaled down by a factor qW < 1. Finally, the improved equation of motion 
(7b) does follow conveniently from a Lagrangian 



!<*>* = OTU 2/ 2 - v W* where V^* = qW V, (7c) 

in sharp contrast to the customary Rayleigh model. Clearly, if V does not contain 
t explicitly the improved energy E\ = mv 2 /2 + V, is approximately conserved, 
unless one goes far from the barrier. 



3.4 The non-Rayleigh branch 

Next, we turn to the very interesting case of the particle climbing up the potential hill in 
figure 1, so that the velocity is essentially antiparallel to the applied force. The particle 
then takes a relatively long time to cover a given distance up the hill. For times large 
compared to I/O,-, employing the asymtotic expansion (A12) we obtain for the velocity 

uM*=-F/mO_ ; ft_fl. (8a) 

Here the dots . . . represent nonleading terms of order F p /mfi^. and the superscript (N)* 
labels the non-Rayleigh branch under consideration. The corresponding acceleration is 
deduced from 



which will not depend on t explicitly if V does not, in sharp contrast to the BCK [Ic, 2a] 
model. The corresponding energy Ej mv 2 /2 + Vj ' becomes an approximate con- 
stant of motion unless we go far from the barrier. 

3.5 Physical interpretation and numerical estimates 

From the conceptual viewpoint we can say that the effect of the yv term in the exact 
frictional equation for dv/dt is roughly simulated by altering the force in the 
dissipationless formulation. The scale factor q^ in eq. (8b) is greater than unity and 
is an increasing function of y. This physically implies that as the coefficient of friction 
increases, the effective potential hill V^- to be surmounted by the particle becomes 
higher and the deceleration grows. Let us estimate a typical numerical range for the ratio 
y/l for which the asymptotic expansion (8a) on the non-Rayleigh branch holds. Suppose 
O -t- 1 w 100, implying that the particle has completely relaxed with respect to the time 
scale fi^ 1 . Relaxation with respect to the other time scale, viz. il~ ! , would also have been 
achieved during this time provided H_? > 1, i.e., 

O_/fi+ > 1/100, i.e. < y/O < 10. (9) 

In other words the analysis based on eqs (8a-c) remains valid over a fairly wide range of 
underdamped as well as overdamped motions, because y/ti = means no-damping and 
10 implies very heavy damping. 



3.6 Comparison with statistical large-friction expansion 

When the input force includes a random component we must examine any possible link 
between our theory and the 'large y expansions' in statistical mechanics, especially in the 
context of the Fokker-Planck [4a, b] equation. It is well known [4a] that the kramers 
distribution function P(x,v,t) in phase space can be systematically represented by a 
series in inverse powers of y (for very large values of y) as 

p = p() + r - 1 P (1) + 0(r- 2 ). (10) 

The first two terms of this expansion, after integration over the velocity, lead directly to 
the Smoluchowski distribution function W(x, t} in configuration space. Reverting to the 
trajectory language one may say that the asymptotic reduction from Kramers to 
Smoluchowski distribution is equivalent to the asymptotic reduction from the underlying 
Langevin [4b] (cf. eqs ((2a), (b)) to the overdamped Rayleigh [4b] (cf. eqs (3b), 7(a)) 
stochastic differential equation (provided the noise is white Gaussian). In this context, the 
main results of our theory (cf. eqs (8), (9)) are entirely new because they pertain to the 
non-Rayleigh branch of motion, the large parameters involved are Q, and the domain of 
validity extends even to the extreme underdamped limit. 
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All our important results are collected in table 6 which should be compared against 
the entries of tables 1 and 2. Some additional points which deserve mention are as 
follows: 

(a) Limited equivalence: We do not claim perfect equivalence between the actual 
frictional eq. (2a) and its proposed dissipationless versions viz. eqs (7b, 8b). As a matter 
of fact, for time-independent static potentials the exact x trajectories would eventually 
come to rest at a potential minimum whereas our x^ trajectories would go on oscillating 
between static confining walls. However, our approximation (8b) is most effective 
when the particle is climbing up the hill in figure 1, and hence it can be profitably 
employed to calculate the first transit time from the bottom to the top of the barrier. Far 
away from the barrier, i.e., in the region of the confining wells, the one-term approxima- 
tions (7b, 8b) may become poor. In these regions nontrivial corrections due to dissipa- 
tion may be systematically made by including the Fp terms of the asymptotic expansions 
(All, 12). 

(b) Time-dependent situations: Of greater physical interest is the case when the 
perturbation force F p contains explicit time-dependence due to applied modulation, 
white noise, etc. It is known that at times large compared to O"" 1 the frictional energy loss 
suffered by the actual jc trajectories tends to be compensated by the energy gain from the 
time-dependent perturbation so that the particle tends to follow the profile of F p (see 
Appendix, eq. (A8)). Obviously, in the absence of net loss our dissipationless formulation 
based on the x^ trajectories should also become equally valid. 



Table 3. Salient features of the dissipationless Newtonian_ model (cf. eqs (7), (8)). 
Useful abbreviations are Q = (Q 2 + }' 2 / 4 )' /2 > ^+ = ft + y/2. fL = Q - y/2, 
qW = ft_/n+ and qW = n+/fi_. 

Item Dissipationless Newtonian model 

Motivation The effect of deleting the yv term in the exact equation of motion is 

sought to be partially restored by suitably scaling the external force 

Parameters (/?)* branch: fi+large, i.e. Q + t > 1 (N)* branch: Q_ large, i.e 

fH_/3>l. Even the limit of zero friction is included in these 
inequalities 

Estimated velocity v^' = F/mSl + , iW = -F/mfi_ 

Effective acceleration d?; w * /At = qW F, duW /d/ = q Wp 

Branch of motion On the improved Rayleigh branch v and / are essentially parallel as 

happens while sliding down the barrier. On the non-Rayleigh branch 

reverse is the case 

Lagrangian (L) w * = mv 2 /2 - qW V, (L) W * = mv 2 /2 - ^' v )' V 

Mechanical energy (){* } * - mv 2 /2 + qW V f , (E) < f Y =mv 2 /2 + q^ t V, which are 

approximately conserved in regions not too far from the barrier 

(i.e., not inside the confining wells) 

The superscript (#)* labels the improved Rayleigh branch and (N)* labels the non-Rayleigh 



quantum wave packet near the barrier, our dissipationless Lagrangian L^ (see eq. (8c)) 
may have bright applicational prospects. Indeed, if the input V represents time- 
independent mean field, the path integral corresponding to the scaled-potential VW can 
be readily computed via the short-time propagator method [5b] of Sethia et al. Next, if V 
includes time-dependent random components, the ensuing dynamics based on yW can 
be treated using the recently developed concept [6] of quantum dissipationless random 
motion. Finally, a detailed application of these ideas to study the quantum stochastic 
resonance phenomenon [7] in a double well potential under the simultaneous influence of 
thermal temperature, linear friction, external noise and sinusoidal modulation is in 
progress and the results will be reported elsewhere. 

Appendix. 

Solution of eq. (6) 

(i) Green 's function. Given the initial condition XQ and VQ we wish to solve eq. (6) written 
in the form 



where l, 0,+ and O_ are given by eqs (5). The differential operator within square 
brackets has the causal Green's function 

-~ -<M'-0 _ e^'-')], (A2 ) 



" ' 20, L 

where 6(t r 7 ) is the unit step function. 

(ii) Formal solution. If the perturbative force F p were absent, the trajectory XB in 
presence of the barrier alone would have been 

I 

XJB=^-[&~ +t (0,-XQ VQ) + e ~ r (fi + ;ico -f- VQ)]. (A3) 

When F p is also present the use of eqs (A2, 3) permits us to convert eq. (Al) into the 
integral equation 

/oo 

x = x B + I dt'G(t, f / )F' /m 
Jo 

\ ( r l , \ r /' 00 1 ^ 

v, */ 1/ L" \J *J t J J 

(A4) 

where F' p = F^xtf), t'} is the perturbative force at the integration time t' when the 
position becomes xtt 1 ). 



Ap = e- f ./ df'e-^F p . (A5) 

Jt 

Then the formal solution (A4) may be recast into the compact form 

(A6) 



20 l|_ m J L m 

(iii) Asymptotic condition. We shall suppose that the static part of the applied force is 
confining at large distances as happens, for instance, in the case of a double-well 
potential. Then, as t oo the position x can not blow up exponentially, implying that the 
coefficient of e n -' in eq. (A6) should vanish, i.e., 



N = 0; VQ l + XQ + VQ- pp/m. (A7) 

(iv) Final expressions. Inserting the conditions (A7) into eq. (A6) we get the following 
exact, general expression for the position: 

x = xoe-^ 1 - [ p - p p e- n+f + A p ]/2mf2. (A8) 

Differentiation with respect to t yields two mutually equivalent representations for the 
velocity, viz., 



v = tt-x - 2f e~ +f + [ p - p p e- +t ]/m (A9a) 

Ap/m. (A9b) 



It is not difficult to verify explicitly that this v satisfies the standard frictional equation 
(2a) of the text along with the specified initial conditions (2b). 

(v) Existence. At this stage some comments on the mathematical validity of the above 
analysis are in order. Although the adjective 'perturbative' was used to describe F' our 
eqs (A8), (9) are exact because no assumption apart from confinement has been made as 
regards the shape and strength of F?. Also, the functions p , A, and A p in eq. (A5) exist 
throughout the interval < t < oo provided the integrand e~~ -''F' remains bounded at 
all finite t 1 and drops faster than 1/V as t 1 -* oo. Finally, although the top of the barrier in 
figure 1 corresponds to a point of unstable equilibrium, our solutions (A8, 9) are quite 
stable with respect to choice of the initial location XQ because of the important identity 
(A7) satisfied by VQ. 

(vi) Lowest order estimates. We shall now evaluate the integrals appearing in eq. (A5) 
approximately without employing any specific model for F f , but making a few reasonable 
assumptions. Suppose that F P (JE, t) is a slowly varying function of space and time so that 
its derivatives may be neglected - a situation which is somewhat reminiscent of the 
potential smoothness assumption in the quantum WKB method. Also, suppose that the 
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particle has relaxed with respect to the time-scale Hi 1 implying that O f 1. Then, 
partial integration in eq. (A5) yields 

. P F P /OH ApFpH-Fp/n_, (A10) 

where terms of order Fp/tij. have been dropped. With the help of these estimates we can 
describe motion on two different types of branches as explained below. 

(vii) The improved Rayleigh branch. If the particle starts sliding down the barrier in 
figure 1 its speed may be relatively substantial. Relaxation with respect to the time scale 
Q" 1 happens soon so that one can take Cl+t 1. Neglecting the e" - 1 ' terms in the first 
representation (cf. eq. (A9a)) for the velocity and using the estimate (A10)"one obtains 



_ , 

m mu) + mi q. 

where the superscript (R}* labels the 'improved' Rayleigh branch and F = ml 2 x + F p is 
the full external force. Note that the customary Rayleigh velocity is just F/my instead (cf. 
eq. (3b)). Other dynamical functions of interest, viz., the effective acceleration, 
Lagrangian and conserved energy associated with v^ are reported explicitly in table 3. 

(viii) The non-Rayleigh branch. If the particle starts climbing up the hill in figure 1 its 
speed may get reduced substantially. Since the particle will take a relatively long time to 
cover a given distance, relaxation with respect to the time scale fil 1 may occur so that we 
can take fi, 1 . From the second representation (cf. eq. (A9b)) of the velocity we obtain 
the asymptotic expansion 

(A12) 



where use has been made of the estimate (A 10), and the superscript (Af)* labels the non- 
Rayleigh branch. Note that vW is outside the scope of the conventional Rayleigh 
velocity. The relevant acceleration, Lagrangian and conserved energy are again displayed 
explicitly in table 3. 
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Abstract. Interaction of Bernstein mode wave with ion-acoustic turbulence is treated as the 
plasma maser effect. Ion-acoustic turbulence is considered as resonant mode, propagating 
orthogonal to the test Bernstein mode. The external magnetic field, in the direction of ion-acoustic 
turbulence, is playing key role in transferring energy in upconversion process from resonant low 
frequency wave to the high frequency wave through modulated electric field. It is shown that the 
direct coupling term and the associated reverse absorbtion effect do not contribute individually to 
the growth or damping to the Bernstein mode. Only the polarization coupling term is found to be 
effective in the energy upconversion process as the external magnetic field is providing momentum 
in the direction of propogation of Bernstein mode. The polarization coupling term is identified to be 
the dominant upconversion factor in the plasma maser effect. 

Keywords. Weak turbulence theory; resonant and nonresonant waves; linear response theory; 
growth rate. 
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1. Introduction 

Anomalous radiation phenomenon in the presence of enhanced low frequency fluctuation 
is often observed in laboratory and space plasma. Such radiation phenomena of high 
frequency wave may be interpreted as an energy upconversion process. The possibility of 
growth of high frequency wave from low frequency fluctuation was indicated by Roses 
[1]. The theory of Plasma maser effect as an effective upconversion process, was put 
forward independently by Nambu [2] and Tsytovich [3]. Plasma maser effect is a mode- 
mode coupling process in which turbulent energy from low frequency wave is transferred 
to the high frequency wave. Since most of turbulent energy in astrophysical plasma is 
contained in low frequency fluctuations, many energy upconversion processes in it can be 
interpreted through plasma maser effect [4]. Recent numerical simulation [5] shows that 
the accelerated electrons by the electrostatic lower hybrid wave effectively generates the 
electromagnetic radiation at the electron gyrofrequency; this suggests that the theory of 
plasma maser effect can be effective for the interpretation of many radiation phenomena 



Plasma maser effect occurs when resonant as well as nonresonant waves are present in 
the system. The resonant waves are those for which the Cherenkov resonant condition is 
satisfied: v k uj = 0; whereas the nonresonant waves are those for which both the 
Cherenkov resonance condition and the nonlinear scattering conditions are not satisfied 
i.e. ft K v ^ and (K k} v (1 uS) ^0. Here, u and O are frequencies of the 
resonant and nonresonant waves respectively, and k and K are the corresponding wave 
numbers. According to the recent weak turbulence theory, the standard mode coupling 
process consists of three parts, (a) the three wave resonance, (b) the nonlinear Landau 
resonance and (c) the plasma maser effect. Plasma maser effect co-exists with the reverse 
absorption effect due to slow change in distribution function. It is effective in strongly 
magnetised plasma where a free source of energy for the particles is available. 
For unmagnetised plasma, plasma maser effect vanishes [6]; in fact it exactly cancels the 
reversed absorption effect. In this paper we consider ion-acoustic wave fluctuation be 
present in the system and consider its interaction with high frequency Bernstein mode 
when they are propagating in direction orthogonal to each other. Thus we observe growth 
of nonresonant Bernstein mode through plasma maser effect. In mis orthogonal 
interaction the transfer of wave energy is possible due to the presence of magnetic field. 

2. Formulations 

We consider a homogeneous plasma imbedded in an external uniform magnetic field BQ 
along z axis. The interaction of high frequency Bernstein mode with ion-acoustic turbu- 
lence is studied using the Vlasov Poisson equations 

^_ - H-L(E ^ x ^ JL f(- - \-Q 

~dt^ V "d?~m\^ + ' c Jdv-' e ^ V '^~ ^ 

f e (r t v,t)dv. (2) 

We consider the Maxwellian distribution for both the background and beam electrons 
given by 

\ 3 /2 



' ' I ' "* X r 1 f \2l I /O\ 

+ (Vn VQ) J , (3) 



where T),, T are temperatures of the background and beam plasma, 6 = w^/no ( <C 1) is 
the ratio of the beam electron density n b to the background plasma density no, VQ is the 
beam electron velocity, (2T/m) 1 ' 2 is the thermal velocity (VQ ^>(2T/m} 1 ' 2 ) and the sub- 
scripts _L and || denote components perpendicular and parallel to the external magnetic 
field. 

tn th<=>. linp.ar-rf=tfinnnsp. fhp.nrv of f\ hirhnlftnt nlflsma T71 the imnfirfnrhp.d 



d _ d 
-z: + V'-^- 
dt dr 



Using Fourier transforms 
A(r,u,0 = Y~ 



e \ -y x BQ d 
.m/ c dv 



m 






VtJ)exp{/(- ? 



(6) 



(7) 



the Fourier component of/i e (fc,cj) from (6) will be 



m 



+ 



(8) 



where /0 + represents small imaginary part of the frequency uj, and k is the wave number 
of ion-acoustic wave fields. 

To this quasisteady state, we apply perturbation fj,8Ei, of test high frequency electro- 
static Bernstein mode wave (/z <C e) such that the total perturbed electromagnetic fields 
and the distribution function are given as 

8E = , 



Sf =? 



and 



where 8E[h, A are the modulation electric fields and 5///,, A/ are electron distribution 
function corresponding to modulation fields. 
Rewriting the Vlasov equation (1), 



vx 



5f 9r m ^ ' 

-L&E(rf) (f 
m ' dv 

We collect the terms of order LJL, ye and 



e -. d 

=8E h -zz 
m ov 



_d_ 

dv 



These respectively are 



P8f lh = -E l ~ 
m ov 



m 



m 




ov 
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We obtain 6f h , 8f ih and A/ by integrating along the unperturbed orbits [8], using 
cylindrical co-ordinate system to the velocity space, as follows. 
From (9) we have 



e h (K,n)=^8E h 



d 



a b 







dv 



(12) 



where a Kv/l e . From (10) we have 



Kf f& ^ , 

- <5//, (J?:, O) -I- 



\K-k\dv\\ 



. (13) 



From (11) we obtain 






a.b 



K, 



d k\\ d 



,\ K ~ 
The Poisson equation gives the modulated field as 

- f 

V SE [h = -4irn e e I f ih (K - k}dv 

so that 

6E [h (K-k} = 



fle(k") 



(14) 
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Now we obtain dielectric response function for high frequency Bernstein mode wa> 
from Poisson equation: 



V-8E h = -4<rren e {6f h + Af}dv, 
where the integration is taken over the entire velocity space, which takes the form 



The dispersion relation e/,(AT, 0) can be written as 

(, n) = eo (K, n) + e rf (^, n) + e/ ,( 

where eo(^, fi) is the linear part, and is given by 



_^ 

" UMe OV\_ 

is the direct coupling part given by 



and e p (^, fit) is the polarization coupling given by 



M '---< 
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Intregrating (17) by parts we have, 

, n) = i + T 






v dv^ \aj + kv + iO + 
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m\ 2 
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(20) 



where / a denotes the modified Bessel function. While integrating (18) by parts we have 
assumed 



Since A and C are real and B oc D* (* denotes complex conjugate) we have 



(21) 



Here Im is used to denote imaginary part. 
After partial integration, we obtain 
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d/dv\\f Qe (v\\) 
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where 



in obtaining 



Im 



k\\ 



dv\\ 



we have used the relation 
1 



Im 



= 7ri6(uj kv} 



and 



Next considering the fact that for Bernstein mode, the most dominant contribution to th 
Bessel's sums comes from the terms a x = 1, we estimate 



- JVT 



where A2 is given by 

yoo ^ 

A 2 = / ^i () 4 2?r ~ / (fc . (v i 
7o ()v -i 



The expression for Im e p (Eq. (19)) now becomes 

, ,4 



3. Instabilities 

We calculate the growth rate using the formula 



ft ftdeo/dft 

where ft r denotes the real part of ft. The second part is due to revesre absorption effect, 
which in our case is 



~ Kj. ^ (ft 2 - 

where 



dt " 

After partial integration this contribution due to the reverse absorption effect becomes 
zero. 
The growth rate due to polarization coupling term is found to be 

{(u/k\\ 1>o) 1 
x ' " ' \ i 94. 1 
7. f \^ ) 
v l J 

where 



is the ion-acoustic turbulence wave energy. Here W denotes electron number density and 
r e denotes electron temperature. This growth rate is +ve for VQ > u/k\\, therefore only 
the beam plasma is taking part in this process. 



4. Discussions 

Plasma maser effect which comes from the imaginary part of the nonlinear dispersion 
relation, Ime/^AT, ft), is accompanied by absorption effect d 2 eo(K,l,fQ e (t})/2dfldt 



condition i.e., k\\v\\ ui = U. in the absence 01 magnetic Held or in a closed system where 
external source of momentum and particle are not available there is exact cancellation 
between plasma maser effect and the reverse absorption effect i.e., Im/,(AT,f2) + (1/2) 
d 2 eo(K, l}/dldt 0. However for a open system with energy and particle supply, the 
cancellation does not occur and plasma maser interaction is effective. 

In the present case the contribution to Im e/, due to the direct coupling part is found to 
be identically zero. As the resonant ion-acoustic wave is propagating orthogonally there 
is no direct transfer of momentum to the Bernstein mode from ion-acoustic turbulence 
and the direct coupling term will not contribute in the upconversion process. For the same 
reason the contribution due to reverse absorption effect is again found to be zero in this 
particular problem. However, in the presence of the external magnetic field, the system is 
equivalent to externally driven rotation and we have non-zero contribution from polariza- 
tion coupling part to the imaginary part of nonlinear dispersion relation Ime/,. 

In the expresion for growth rate (eq. 24), we have only surviving +ve contribution due 
to the polarization coupling term for VQ > w/k\\. Here, the ambient magnetic field carries 
momentum in the perpendicular to its direction, i.e., towards the propagation of non- 
resonant Bernstein mode in the system. From (eq. 22), for the external magnetic field 
B which gives $l e > we have Im e p > and the plasma maser effect becomes 
ineffective. This clarifies the role of external magnetic field in this energy upconversion 
process. 

In earlier studies [9] on plasma maser effect where interaction of two electrostatic 
waves is considered, both the direct coupling part and the polarization coupling part of 
the nonlinear dispersion relation contributed to the growth rate. In that, the interacting 
waves viz. the resonant ion-acoustic turbulence wave and the nonresonant Langmuir 
wave had common component of their propagation vector in the z-axis so that there is 
direct exchange of momentum among these waves. 

The electron cyclotron maser mechanism [10, 11] is driven by a loss cone type of 
anisotropy. The cyclotron maser mechanism would work if the relativistic electron satisfy 
fiAT[|-y|| - fi e /7 = where ATy is parallel component of high frequency electromagnetic 
radiation and 7 = (1 v 2 /c 2 }~~ ] ' . The relativistic correction term provides a mechanism 
which couples the perpendicular motion v] with waves and allows the free energy perpen- 
dicular to the magnetic field to be transferred to electromagnetic radiation. Again the 
standard cyclotron maser mechanism predicts that the radiation mechanism is effective 
only at frequencies just above the cutoff frequency of high frequency electromagnetic 
wave which is very close to electron cyclotron frequency l e . However the nonrelativistic 
theory of cyclotron maser theory proposed by Sagdeev and Sharfranov [12] states that 
cylotron maser instability can develop in nonrelativistic magnetised plasma with aniso- 
tropic distribution function. The cyclotron maser instability cannot developie [13] in 
tenuous plasma N w l,(Nis refractive index) for the longitudinal propagation. Further it 
has been established with the relativistic consideration [14] that cyclotron maser instabi- 
lity can developie at N > 1 if only the distribution function is anisotropic. The instability 
disappears when W 1, and the negative aborption is exponentially small at N 2 1 <C 1. 

On the other hand in plasma maser effect, the free energy to be transferred to the high 



nonrelativistic electrons. It does not require any frequency matching condition among me 
interacting waves; so for the case of electromagnetic radition [15] in the presence of low 
frequency turbulence field the plasma maser is effective well above the local electron 
gyrofrequency. On the other hand the standard cyclotron maser mechanism [10,11] 
predicts that the generation is effective only at frequencies just above the cutoff frequency 
of electromagnetic radiation which is very close to JT2 tf . In plasma maser effect the 
accompanying reverse absorption effect competes with the direct coupling part of the 
dispersion relation. The dominant part in plasma maser effect comes from polarization 
coupling part. The plasma maser instability becomes ineffective as the strength of the 
external magnetic field decreases. 
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